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FOREWORD 


Atmospheric  parameters  are  essential  to  the  research  and  development  of 
missiles  and  aerospace  vehicles.  In  the  early  1960's,  the  need  was  recognized 
for  realistic  atmospheric  models  derived  in  a  consistent  manner  for  each  of  the 
several  major  test  ranges.  An  atmospheric  model  derived  from  statistical  data 
for  a  particular  geographical  location  is  referred  to  as  a  reference  atmosphere. 

The  first  Range  Reference  Atmosphere  (RRA)  was  issued  in  1 §63  by  the 
Inter-Range  Instrumentation  Group  (IRIG)  for  Cape  Kennedy,  Florida,  and  was 
followed  by  additional  publications  for  several  ranges  up  to  1974.  Since 
that  time,  improved  upper  air  data  bases  have  become  available  from  which  to 
develop  the  RRA.  These  resulted  from  the  extended  period  of  records  and  from 
improvement  in  the  upper  air  measuring  program  by  rocketsondes  for  altitudes 
above  the  rawinsonde  ceiling  of  30  km.  Revised  and  improved  RRAs  are  justi¬ 
fied  for  the  following  reasons: 

1)  Needs  for  more  definitive  statistical  atmospheric  models  have  arisen 
because  of  changes  and  advances  in  aerospace  technology.  The  Space  Transpor¬ 
tation  System  (Space  Shuttle)  is  one  example. 

2)  Most  ranges  now  have  an  extended  and  improved  upper  air  data  base 
from  which  to  develop  a  more  definitive  RRA. 

3)  There  are  requirements  for  RRAs  for  new  ranges  and  range  sites. 

4)  There  have  been  scientific  advances  in  understanding  the  upper 
atmospheric  structure  and  physical  relationships. 

5)  Advances  in  statistical  modeling  techniques  have  been  made  because 
of  the  general  availability  of  high-speed  electronic  computers.  These  have 
led  to  the  adoption  of  advanced  concepts  in  atmospheric  modeling. 

For  these  reasons,  the  Range  Reference  Atmosphere  Committee  (RRAC)  was 
tasked  by  the  Range  Commanders  Council  Meteorology  Group  (RCC  MG)  to  estab¬ 
lish  new  and  improved  RRAs.  The  purpose,  scope,  and  objectives  of  this  task 
are  outlined  in  the  following  paragraphs. 

Purpose:  This  committee.  Task  MG-1 ,  establishes  RRAs  for  the  several 
ranges  as  provided  by  the  RCC.  An  RRA  is  a  model  of  the  Earth's  atmosphere 
over  a  geographical  location  of  interest,  for  use  by  DOD  and  other  U.S. 
Government  range  users.  The  RRA  is  used  to  provide  planning  data  for  eval¬ 
uating  environmental  constraints  for  the  particular  configurations  of 
environment-sensitive  systems  and  components  being  developed  or  undergoing 
tests. 

Scope:  Using  the  best  available  upper  atmosphere  data  base  to  include 
rawinsonde,  rocketsonde  and  possibly  other  high-altitude  data  sources  for 
the  range  location,  the  task  is  to  establish  a  model  of  certain  statistics 
for  wind  and  thermodynamic  quantities  derived  in  a  uniform  manner  and  pub¬ 
lished  in  a  standardized  format. 


N'.-. 


vi  i 


Objectives:  The  wind  statistics  shall  be,  insofar  as  practical,  modeled 
to  be  consistent  with  rigorous  mathematical  probability  properties  of  the  mul¬ 
tivariate  normal  probability  theory.  The  thermodynamic  quantities  statistics 
shall  be,  insofar  as  practical,  modeled  to  be  consistent  with  the  hydrostatic 
equation,  the  equation  of  state,  and  the  probability  principles  that  are 
related  through  these  physical  equations.  The  document  shall  serve  as  an 
authoritative  source  of  information  and  as  an  atmospheric  model  for  a 
particular  range.  The  first  in  the  series  of  revised  RRAs  to  be  published 
is  for  Kwajalein  Missile  Range  (KMR)  (publication  date  December  1982).  The 
altitude  range  required  for  KMR  is  0  to  70  km.  The  order  of  priority  for 
the  subsequent  publications  is: 

Range  Altitude  Range  Required 


1. 

AFFTC/Edwards  AFB,  CA 

0  -  70  kma 

2. 

ESMC/Cape  Canaveral  AFS,  FL 

0  -  70  km 

3. 

WSMC/Vandenberg  AFB,  CA 

0  -  70  kma 

4. 

WSMR/Uhite  Sands,  NM 

0  -  70  km 

5. 

PMTC/Point  Mugu,  CA 

0  -  70  km 

6. 

UTTR/Dugway  (Michael  AAF),  UT 

0  -  30  km* 

7. 

AD/Eglin  AFB,  FL 

0  -  30  km 

8. 

ESMC/Ascenslon  Island 

0  -  70  km  (Terminates  at  66  km 
because  of  insufficient  data) 

9. 

NASA/Wallops  Flight  Center,  VA 

0  -  70  km 

10. 

Taquac  (Guam) 

0  -  30  km 

11. 

PMTC/Barking  Sands,  HI 

0  -  70  km 

In  keeping  with  the  RCC's  objective  of  standardization,  the  modeling 
techniques,  basic  text,  and  tabulation  format  are  to  be  the  same  for  all  RRAs. 
These  new  and  revised  RRAs  present  not  only  the  mean  values  of  the  thermody¬ 
namic  quantities  (pressure,  temperature,  virtual  temperature,  and  density), 
but  also  Include  statistical  measures  for  the  dispersion  (l.e.,  standard 
deviations  and  skewness  coefficients).  New  quantities  presented  are  water 
vapor  pressure  and  dewpoint  temperature.  The  statistical  modeling  for  the 
wind  Is  entirely  new.  The  new  approach  uses  the  properties  of  the  bivariate 
normal  probability  distribution  function. 


a.  Use  rooketsoncf  data  from  PMTC/Point  Mugu  for  altitudes  above  30  km. 

b.  Consider  augmenting  data  base  from  Ely  or  Salt  Lake  City. 
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All  final  computations  were  performed  by  the  United  States  Air  Force 
Environmental  Technical  Applications  Center  (USAFETAC)  in  response  to  a  task 
from  Eastern  Space  and  Missile  Center  (ESMC). 

The  text  was  prepared  jointly  by  USAFETAC  and  the  NASA/George  C.  Marshall 
Space  Flight  Center's  Space  Sciences  Laboratory,  Atmospheric  Sciences  Division. 
The  editing  and  preparation  of  the  draft  manuscript  were  performed  by  the 
NASA/MSFC  organization. 

The  cochairmen  express  their  gratitude  to  all  RRAC  members  and  their 
respective  colleagues  who  have  made  significant  technical  contributions  to 
the  establishment  of  these  RRAs. 

Special  thanks  are  tendered  to  Lt.  B.  Novograd  for  his  dilligence  in 
forming  the  many  computations  and  the  development  of  the  primary  tables,  I 
through  IV.  Special  thanks  goes  to  Lt.  F.  Wirsing  for  editing  and  formulating 
the  equations  for  the  derivable  thermodynamic  equations.  These  gentlemen  per¬ 
formed  this  outstanding  work  under  the  direction  of  Major  B.  Lilius,  USAFETAC. 

Grateful  acknowledgment  goes  to  Mrs.  Annette  Tingle,  NASA/MSFC,  for  edit¬ 
ing  the  draft  manuscript. 

The  RRAC  consists  of  representatives  from  the  U.S.  Air  Force,  U.S.  Army, 
National  Aeronautics  and  Space  Administration,  U.S.  Navy,  and  National  Oceanic 
and  Atmospheric  Administration.  The  committee  members  for  the  RRA  for  the 
first  publication  are: 

G.  G.  Boire,  WSMC 

0.  H.  Daniel,  ESMC 

R.  de  Viol  ini,  PMTC 

F.  G.  Finger,  NOAA/NWS 

E.  E.  Fisher,  HQ  AFSC 

B.  R.  Hixon,  PMTC 

J.  M.  Hobbie,  KMR 

E.  J.  Keppel ,  AD 

S.  F.  Kubinski,  WSMR 

F.  J.  Schmidlin,  NASA/WFC 

0.  E.  Smith  Maj.  B.  W.  Galusha 

Cochairman,  NASA/MSFC  Cochairman,  USAF/ETAC 
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CHAPTER  I.  INTRODUCTION 


A.  Definition  and  Purpose  of  the  Range  Reference  Atmosphere 
A.l  Definition 

A  reference  atmosphere  is  a  statistical  model  of  the  Earth's  atmosphere 
derived  from  upper  air  measurements  over  a  particular  geographical  location. 
Hence,  these  Range  Reference  Atmospheres  (RRAs)  are  atmospheric  models  devel¬ 
oped  by  the  Range  Reference  Atmosphere  Committee  (RRAC)  in  response  to  a  task 
by  the  Range  Commanders  Council  Meteorology  Group  (RCC  MG)  and  published  by  the 
RCC  Secretariat.  The  RCC  MG,  formerly  called  the  Inter-Range  Instrumentation 
Group/Meteorology  Working  Group  (IRIG/MWG),  published  a  series  of  RRAs  during 
the  period  1963  through  1974. 

A.  2  Purpose 

A  series  of  revised  and  expanded  RRAs  are  to  be  published  for  locations 
of  interest  to  the  RCC.  These  publications  are  to  serve  as  authoritative 
reference  sources  on  certain  upper  air  statistics  and  as  atmospheric  models 
for  particular  range  sites.  The  technical  usefulness  of  these  documents  for 
the  ranges,  range  users,  U.S.  aerospace  industries,  and  the  scientific  commu¬ 
nity  is  recognized  because  of  the  standardization  of  the  development  tech¬ 
niques  and  the  presentation  of  the  tabulations. 

B.  Scope  of  the  Range  Reference  Atmosphere  and  Arrangement  of  Tables 
B.T  Scope 

^The  RRA  contains  tabulations  for  monthly  and  annual  means,  standard  devia¬ 
tions,  and  skewness  coefficients  for  windspeed,  pressure,  temperature,  density, 
water  vapor  pressure,  virtual  temperature,  and  dewpoint  temperature;  the  means 
and  standard  deviations  for  the  zonal  (U)  and  meridional  (V)  wind  components; 
and  the  linear  (product  moment)  correlation  coefficient  between  the  wind  com¬ 
ponents.  These  statistical  parameters  are  tabulated  at  the  station  elevation, 
at  1-km  intervals  from  sea  level  to  30  km,  and  at  2-km  intervals  from  30  to 
90  km.  The  wind  statistics  are  given  at  approximately  10  m  above  the  station 
elevations  and  at  altitudes  with  respect  to  mean  sea  level  thereafter.  For 
those  range  sites  without  rocketsonde  measurements,  the  RRAs  terminate  at 
30  km  altitude,  or  they  are  extended,  if  required,  when  rocketsonde  data 
from  a  nearby  launch  site  are  available.  There  are  four  sets  of  tables  for 
each  of  the  12  monthly  reference  periods  and  the  annual  reference  period. 

B.2  Arrangement  of  Tables 

The  statistical  parameters  for  the  RRA  models  are  presented  in  four  tables, 
as  outlined  in  the  following  paragraphs. 

Table  I  contains  all  the  wind  statistical  parameters.  This  table  gives 
the  monthly  and  annual  means  and  standard  deviations  of  the  U  and  V  wind  com¬ 
ponents  and  the  linear  (product  moment)  correlation  coefficient  between  these 


two  components;  the  mean,  standard  deviation  and  skewness  coefficient  of  the 
windspeed;  and  the  number  of  wind  observations  (sample  size). 

Table  II  contains  the  monthly  and  annual  means,  standard  deviations,  and 
skewness  values  of  pressure,  temperature,  and  density,  and  the  number  of  obser¬ 
vations  used  for  each  of  these  thermodynamic  quantities. 

Table  III  contains  the  monthly  and  annual  means,  standard  deviations  and 
skewness  values  of  the  water  vapor  pressure,  virtual  temperature  and  dewpoint, 
and  the  number  of  observations  for  each  of  these  moisture-related  quantities. 

The  statistical  parameters  for  water  vapor  pressure  and  dewpoint  terminate  at 
15  km  altitude.  Above  15  km  the  statistical  parameters  for  virtual  temperature 
are  considered  to  be  the  same  as  those  for  temperature. 

Table  IV  contains  the  monthly  and  annual  mean  atmospheric  models  for  the 
thermodynamic  variables:  pressure,  virtual  temperature,  and  density.  This 
table  is  derived  from  the  monthly  and  annual  mean  virtual  temperature  versus 
altitude  (geometric)  using  the  hydrostatic  equation  and  the  equation  of  state. 
Also  presented  is  the  geopotential  height  corresponding  to  the  tabulated  geo¬ 
metric  altitudes. 

The  physical  unit  for  all  wind  parameters  is  meters  per  second.  The  phy¬ 
sical  unit  for  pressure  is  millibars;  for  temperature  and  virtual  temperature, 
degrees  Kelvin;  for  density,  grams  per  cubic  meter;  and  for  water  vapor  pressure, 
millibars.  In  all  cases  the  skewness  coefficient  and  the  correlation  coefficient 
between  wind  components  are  unitless.  All  reference  to  altitude  is  geometric 
altitude  and  is  expressed  in  kilometers.  All  reference  to  height  is  geopoten¬ 
tial  height  and  has  the  unit  geopotential  meters  or  kilometers.  All  geometric 
altitudes  and  geopotential  heights  are  with  respect  to  mean  sea  level. 

C.  Data  Quality  Control  Procedures 

A  small  portion  (less  than  10  percent)  of  the  soundings  in  the  data  base 
used  to  cal-uiate  the  RRA  tables  contained  erroneous  data  values.  The  soundings 
which  contained  these  erroneous  values  were  eliminated  from  the  data  base  using 
the  following  procedures: 

1)  Soundings  containing  gaps  in  their  height  data  greater  than  200  mb  were 
rejected.  This  step  was  taken  because  some  soundings  only  contained  height 
values  at  their  "mandatory"  pressure  levels,  which  were  occasionally  missing, 
resulting  in  soundings  with  no  height  information  at  all. 

2)  An  initial  set  of  RRA  statistics  was  computed  using  all  the  remaining 
soundings.  This  initial  set  of  statistics  was  used  to  determine  data  limits 
for  the  temperature,  pressure,  U  and  V  components  of  the  wind,  and  the  dewpoint 
(for  the  0-  to  30-km  portion  of  the  RRA)  or  the  density  (for  the  30-  to  90-km 
portion  of  the  RRA).  The  lower  (upper)  data  limits  were  set  at  the  mean  value 
for  a  specific  parameter,  minus  (plus)  six  standard  deviations  of  that  quantity. 
One  pair  of  data  limits  was  computed  for  each  of  these  parameters:  month  of  the 
year  and  data  level. 
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3)  This  initial  set  of  data  limits  was  then  used  to  screen  the  data  base. 
All  the  soundings  that  contained  values  outside  these  data  limits  were  rejected. 
A  new  RRA  was  then  computed  using  the  screened  data  base.  This  second  RRA  was 
used  to  generate  a  second  set  of  data  limits. 

4)  The  second  set  of  data  limits  was  then  used  to  screen  the  data  base 
further.  A  new  RRA  was  again  generated.  The  skewness  values  in  this  RRA  were 
then  evaluated,  according  to  empirical  criteria  specified  in  section  II. A. 3  of 
this  document  for  the  winds,  and  according  to  criteria  in  section  III. A. 3  for 
the  thermodynamic  quantities.  If  these  criteria  were  satisfied,  the  new  RRA 
was  then  used  to  generate  a  final  set  of  data  limits,  which  were  used  to 
control  the  quality  of  the  data  base  for  the  final  version  of  the  RRA. 

5)  Occasionally,  the  third  RRA  that  was  generated  did  not  satisfy  all 
of  the  skewness  criteria.  This  indicated  that  some  incorrect  values  were 
still  present  in  the  data  base.  To  complete  quality  control,  steps  3  and  4 
were  repeated  for  additional  iterations  (usually  one  or  two)  until  the  result¬ 
ing  RRA  satisfied  the  skewness  criteria.  At  that  point,  a  final  set  of  data 
limits  was  generated.  This  final  set  of  data  limits  was  then  used  to  control 
the  quality  of  the  data  base  and  generate  the  final  RRA. 

D.  Organization  of  the  Chapters 

Because  there  are  plans  to  publish  a  series  of  RRAs,  comments  on  the 
special  organization  of  the  document  are  in  order.  The  RRA  document  is 
arranged  in  four  chapters.  Chapter  I  is  the  introduction.  Chapter  II, 

Wind  Statistics  and  Models,  contains  the  techniques  used  to  arrive  at  the 
wind  statistical  parameters,  table  I,  and  the  probability  functions  that 
are  to  be  used  as  wind  models  to  derive  several  wind  statistics.  Chapter 
III,  Statistics  of  Thermodynamic  Quantities  and  Models,  contains  the  tech¬ 
niques  used  to  arrive  at  the  thermodynamic  and  moisture- related  statistical 
parameters  given  in  tables  II  and  III  and  the  atmospheric  thermodynamic 
model  presented  in  table  IV.  This  chapter  also  contains  sets  of  equations 
to  calculate  several  atmospheric  properties.  Chapter  IV  contains  the  gen¬ 
eral  conclusions  and  recommendations.  These  four  chapters  are  reprinted 
without  change  for  each  documented  RRA  to  assure  consistency  and  for  expe¬ 
diency  in  preparing  the  documentation.  To  account  for  variations  particular 
to  a  specific  RRA,  two  appendixes  have  been  included.  Appendix  A,  Examples 
of  Wind  Statistics,  is  designed  to  give  a  few  illustrative  examples  of  wind 
statistics  for  the  specific  RRA  and  cursory  observations,  comparisons,  or 
comments  on  wind  statistics.  Appendix  B,  Range  Specific  Information,  is 
designed  to  present  specific  information  particular  to  the  range,  such  as 
geographical  location,  data  base,  etc.,  and  any  cursory  observations  or 
comments  on  the  thermodynamic  quantities. 

Read  these  appendixes!  They  are  located  as  the  last  two  units  in  the 
document  because  they  may  vary  in  length  depending  on  the  circumstances. 
Appendixes  A  and  B  and  tables  I,  II,  III,  and  IV  are  the  only  differences 
among  the  RRA  documents  published  in  this  new  RRA  series. 


CHAPTER  II.  WIND  STATISTICS  AND  MODELS 


A.  General  Considerations 
A.l.  Objectives 

An  objective  of  the  RRA  is  to  furnish  minimum  tabulation  for  the  wind 
statistics.  To  meet  this  objective,  the  bivariate  normal  probability  distri¬ 
bution  was  adopted  as  a  statistical  model  for  the  wind  treated  as  a  vector 
quantity  at  the  RRA  data  levels.  Only  five  statistical  parameters  are  required 
to  completely  describe  this  probability  function.  In  Cartesian  coordinates 
these  parameters  are  the  means  and  standard  deviations  of  the  two  orthogonal 
components  and  the  correlation  coefficient  between  the  two  components.  These 
five  statistical  parameters  for  the  U  and  V  (meteorological  coordinates)  com¬ 
ponents  are  given  in  table  I.  The  statistical  properties  of  the  bivariate 
normal  probability  distribution  are  used  to  derive  many  wind  statistics  that 
are  of  interest  to  the  ranges  and  range  users.  This  procedure  produces  con¬ 
sistent  wind  statistics  that  are  connected  through  rigorous  mathematical 
probability  functions.  By  using  these  functions,  extensive  tabulations  of 
wind  statistics  are  avoided. 

The  statistical  properties  of  the  bivariate  normal  probability  distribu¬ 
tion  presented  for  the  vector  wind  statistical  model  are: 

1)  The  wind  components  are  univariate  normally  distributed. 

2)  The  conditional  distribution  of  one  component  given  a  value  of  the 
other  component  is  univariate  normally  distributed. 

3)  The  windspeed  is  of  the  form  of  a  generalized  Rayleigh  distribution. 

4)  The  frequency  distribution  of  wind  direction  can  be  derived. 

5)  The  conditional  distribution  of  windspeed  given  a  value  of  wind 
direction  (wind  rose)  can  be  derived. 

6)  The  five  tabulated  wind  statistical  parameters  with  respect  to  the 
meteorological  U  and  V  coordinate  system  can  be  derived  for  any  arbitrary 
rotation  of  the  orthogonal  axes. 

The  probability  distribution  functions  and  sets  of  equations  to  derive 
wind  statistics  for  the  previously  stated  properties  of  the  vector  wind  model 
are  presented  in  this  chapter.  Symbols  used  are  summarized  in  table  A. 
Illustrative  examples  are  presented  in  appendix  A.  No  attempt  is  made  to 
give  the  derivation  of  the  probability  functions.  The  reader  is  referred 
to  Smith  (1976)  for  some  derivations  and  several  applications  of  the  prob¬ 
ability  distribution  properties  for  wind  statistics. 

A. 2.  Data  Quality  Control 

The  U  and  V  components  of  the  wind  were  used  to  generate  data  limits  set 
at  plus  and  minus  six  standard  deviations  from  the  mean  for  each  of  the 


TABLE  A.  LIST  OF  SYMBOLS  USED  IN  CHAPTER  II 


N  -  The  number  of  wind  measurements  in  table  I 

r  -  A  general  variable  for  the  bivariate  normal  probability  distri¬ 
bution  in  polar  coordinates 

R  -  A  generalized  Rayleigh  variable  used  for  derived  windspeed 
probability  distribution 

R  (U,  V)  -  The  linear  (product  moment)  correlation  coefficient 
between  the  zonal  and  meridional  wind  components  in  table  I 

SK  (W)  -  Skewness  parameter  for  windspeed  in  table  I 

S  (U)  -  The  standard  deviation  of  the  zonal  wind  component  in 
table  I 

S  (V)  -  The  standard  deviation  of  the  meridional  wind  component 
in  table  I 

S  (W)  -  The  standard  deviation  of  windspeed  in  table  I 
t  -  A  standardized  normal  variate  used  in  text  table  B 
U  -  The  zonal  wind  component 

UBAR  -  The  mean  value  of  the  zonal  wind  component  in  table  I 

Y  -  The  meridional  wind  component 

VBAR  -  The  mean  value  of  the  meridional  wind  component  in 
table  1 

W  -  Windspeed  or  modulus  of  wind  vector,  a  scalar  quantity 
WBAR  -  The  mean  value  of  windspeed  in  table  I 
X  -  A  general  component  variable  or  coordinate  axis 

Y  -  A  general  component  variable  or  coordinate  axis 

X  -  A  general  component  mean  value  in  the  [x,y]  coordinate  system 

Y  -  A  general  component  mean  value  in  the  [x,y]  coordinate  system 

a  (alpha)  -  Rotation  angle  for  the  [x,y)  coordinate  system 


TABLE  A.  (concluded) 


0  (theta)  -  Wind  direction  in  the  polar  coordinate  system 

j  (Lambda)  -  A  parameter  in  the  bivariate  normal  probability 
distribution  in  text  table  C 

£  (Xi)  -  The  mean  value  in  the  standardized  normal  probability 
distribution  used  in  text  table  B 

w  (Pi)  -  Constant  =  3.14159  ... 

p  (Rho)  -  The  general  linear  correlation  coefficient  between  the 
two  component  variables  in  the  [x,y]  coordinate  system 

°x>ay  ~  The  general  standard  deviations  of  the  x  and  y  compo¬ 
nent  variables  in  the  [x,y]  coordinate  system. 


quantities.  These  data  limits  were  used  to  screen  the  wind  data  base,  as 
described  in  section  I.C.  The  data  base  was  considered  to  be  free  from  errors 
under  the  following  conditions: 

1)  The  skewness  of  the  windspeed  was  below  4.0  at  data  levels  where  the 
mean  windspeed  was  less  than  15  m/s,  and 

2)  The  skewness  of  the  windspeed  was  below  2.5  at  data  levels  where  the 
mean  windspeed  was  greater  than  15  m/s. 

A.  3  Limitations 

For  the  wind  statistics,  the  correlation  coefficients  for  like  wind  com¬ 
ponents  and  unlike  wind  components  between  altitude  levels  were  not  computed. 
Therefore,  wind  statistics  with  respect  to  altitude  (profile)  cannot  be 
derived  from  the  RRA  statistics.  For  wind  profile  modeling  techniques  the 
user  is  referred  to  Smith  (1976).  However,  the  wind  statistics  at  discrete 
altitudes  are  valid;  all  of  the  probability  distribution  functions  given  in 
chapter  II  can  be  derived  from  the  five  wind  component  statistical  parameters 
contained  in  table  I,  and  the  derived  distributions  can  be  considered  as  wind 
models  at  discrete  altitudes. 

By  convention,  in  the  statistical  literature  Greek  letters  are  used  for 
population  or  theoretically  known  parameters,  and  sample  estimates  are  denoted 
by  English  alphabetical  letters  or  with  a  “hat"  (~)  over  the  Greek  letters. 

In  chapter  II  Greek  letters  are  used  for  the  variances  and  the  linear  corre¬ 
lation  coefficient,  and  the  means  are  denoted  by  X  and  Y  when  dealing  with  the 
bivariate  normal  distribution.  It  will  always  be  understood  that  table  I  con¬ 
tains  sample  estimates  of  the  statistical  parameters  and  they  are  with  respect 
to  the  meteorological  U  and  V  coordinate  system. 

B.  Coordinate  System  and  Computation  of  Statistical  Parameters 
B.l.  Coordinate  System 

Wind  measurements  are  recorded  in  terms  of  magnitude  and  direction.  The 
wind  direction  is  measured  in  degrees  clockwise  from  true  north  and  is  the 
direction  from  which  the  wind  Is  blowing.  The  wind  magnitude  (the  modulus  of 
the  vector)  is  the  scalar  quantity  and  is  referred  to  as  windspeed  or  scalar 
wind.  A  statistical  description  that  accounts  for  the  wind  as  a  vector  quan¬ 
tity  is  appropriate  and  requires  a  coordinate  system. 

For  the  RRA  the  standard  meteorological  coordinate  system  has  been  chosen 
for  the  wind  statistics,  all  tables  of  statistical  parameters,  and  related  dis 
cussions  because  the  coordinate  system  used  in  aerospace  and  related  applied 
fields  has  not  always  been  consistent. 

Using  figure  1,  the  polar  and  Cartesian  forms  for  the  meteorological 
coordinate  system  are  defined: 
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W  =  windspeed,  scalar  wind,  or  magnitude  of  the  wind  vector  in  meters  per 
second . 

e  =  wind  direction.  0  is  measured  in  degrees  clockwise  from  true  north 
and  is  the  direction  from  which  the  wind  is  blowing. 

U  =  zonal  wind  component,  positive  west  to  east,  in  meters  per  second. 

V  *  meridional  wind  component,  positive  south  to  north,  in  meters  per 
second. 

The  components  0  and  W  define  the  polar  form,  and  the  U-V  components 
define  the  Cartesian  forms: 

U  =  -W  sin©  ,  0  £  0  1  360°  (1) 

V  =  -W  cos0.  (2) 

north 


Figure  1.  The  meteorological  coordinate  system. 


It  is  helpful  to  note  the  difference  between  the  mathematical  convention 
for  a  vector  direction  and  the  meteorological  convention  for  wind  direction: 

e  met  =  270  -  e  math  (3) 

when  0  ^  o  math  <  270° 

9  met  =  360  +  (270  -  9  math) 

when  270  <  9  math  £  360° 

B.2  Computation  of  Statistical  Parameters 

The  wind  statistical  parameters  in  table  I  for  the  means  and  standard 
deviations  of  the  U  and  V  wind  components  and  windspeed  and  the  skewness 
parameter  of  windspeed  were  computed  using  the  sums  technique  presented  in 
chapter  III.C.3.  In  addition,  the  linear  (product  moment)  correlation  coef¬ 
ficient  between  the  U  and  V  wind  components,  r  (u,v)  In  table  I,  was  computed. 
This  correlation  coefficient  is  defined  as 


r  (u.v) 


(4) 


These  statistical  parameters  are  with  respect  to  the  Standard  Meteorological 
Coordinate  System. 

C.  Statistical  Wind  Models 


C.l.  Wind  Component  Statistics 


The  univariate  normal  (Gaussian)  probability  distribution  function  is  used 
to  obtain  wind  component  statistics.  In  generalized  notations,  this  probability 
density  function  (pdf)  is 


f  (t)  = 


(5) 


where  t  -  X  -  c/a  is  the  standardized  variate,  with  £  defining  the  mean 
and  ox  the  standard  deviation.  The  probability  distribution  function  (PDF)  is 


X 

F( X )  =  /  f(t)  dt 


(6) 
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Because  this  integral  cannot  be  obtained  in  closed  form,  it  is  widely 
tabulated  for  zero  mean  and  unit  standard  deviation.  For  a  convenient 
reference  for  the  RRA,  selected  values  of  F(X)  are  given  in  table  B.  To 

emphasize  the  connotation  of  probability,  F(X)  is  shown  in  table  B  as  pjxj. 

The  t  values  in  table  B  are  used  as  multiplier  factors  to  the  standard  devia¬ 
tion  to  express  the  probability  that  a  normally  distributed  variable,  X,  is 
less  than  or  equal  to  a  given  value  as 


>{x  ^ 


mean  +  t  o  l  = 


l  =  probability ,  p 


(7) 


For  example,  when  t  =  1.6449,  the  probability  that  X  is  less  than  or  equal  to 
the  mean  plus  1.6449  standard  deviations  is  0.95.  That  value  of  X  that  is 
less  than  or  equal  to  the  mean  plus  1.6449  standard  deviations  is  called  the 
95th  percentile  value  of  X.  Also  given  in  table  B  are  the  numerical  values 
to  express  the  probability  that  X  falls  in  the  interval  X^  and  X^;  i.e.. 


pj  Xj  .‘.X  <  X2  |  =  Interpercentile  Range  ,  ^ 

where 

X1  =  *  1  x 

X  X  H  t  Ox 

For  t  =  1.9602  the  probability  that  X  lies  in  the  interval  X^  and  X2  is  0.95. 

The  values  of  X^  and  i*1  this  example  comprise  the  95th  interpercentile  range. 

For  a  normally  distributed  variable,  the  mode  (most  frequent  value)  and  the 
median  (50th  percentile  value)  are  the  same  as  the  mean  value.  The  means  and 
standard  deviations  of  the  U  and  V  wind  components  from  table  1  are  used  in 
equations  (7)  and  (8)  to  compute  the  percentile  values  and  interpercentile 
ranges  of  the  U  and  V  wind  components.  When  equation  (7)  is  illustrated  on 
a  normal  probability  graph,  a  straight  line  is  formed. 

C.2.  The  Vector  Wind  Model 

Because  wind  is  a  vector  quantity  having  direction  and  magnitude  that  can 
be  expressed  as  two  components  in  an  orthogonal  coordinate  system,  a  probability 
model  that  describes  the  joint  relationship  is  the  bivariate  normal  probability 
distribution.  In  general  component  notation,  the  bivariate  normal  probability 
density  function  (BNpdf)  is 


TABLE  B.  VALUES  OF  t  FOR  STANDARDIZED  NORMAL 
(UNIVARIATE)  DISTRIBUTION  FOR  PERCENTILES 
AND  INTERPERCENTILE  RANGES 


f(X,Y)  = 


2TOxOy  /l  ~  p2 


exp 


2(1  -  p‘) 


(X  -  X)2 


_  2p(X  -  X)  (Y  -  Y)  ,  (Y  -  Y)2 


o  o 
x  y 


£  X  <.  »  and 


-  00  £  Y  £  OO  , 


(9) 


where  the  five  parameters  are  x,y,  the  component  means;  ox>  cy  the  component 

standard  deviations;  and  p,  the  correlation  coefficient  between  the  two  compo¬ 
nent  variables,  X  and  Y. 

For  many  applications  the  interest  is  in  determining  the  probability  that 
a  point  |x,y|  will  fall  within  a  contour  of  equal  probability  density-  The 
exponential  terms  of  equation  (9),  when  set  equal  to  a  constant,  A2,  give  a 
family  of  ellipses  depending  on  the  value  of  the  constant.  The  ellipses  have 

a  conroon  center  at  the  point  {xj}.  Integration  of  equation  (9)  over  the 
region  bounded  by  the  contours  of  equal  probability  density  gives 

-X2 

P( A)  =  1  -  e  2(1  '  (,2)  .  (10) 

Solving  for  A2  and  replacing  P(a)  by  p  gives 

A2  =  -  2  (1  -  p2)  In  (1  -  P)  .  (11) 

Now  define 

Ae  =  /T  /-  In  (1  -  p)  .  (12) 

For  ready  reference  and  comparisons,  Ag  is  shown  in  table  C  for  selected 
values  of  p. 
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TABLE  C.  VALUES  OF  X  FOR  BIVARIATE  NORMAL 
DISBRIBUTION  ELLIPSES  AND  CIRCLES 


p(%) 

X 

e 

(ellipse) 

X 

c 

( circle) 

P(%) 

X 

e 

(ellispe) 

X 

c 

(circle) 

0.000 

0. 0000 

0.0000 

65.000 

1.4490 

1.0246 

5.000 

0. 3203 

0.2265 

68.268 

1. 5151 

1.0713 

10.000 

0.4590 

0. 324G 

70.000 

1.5518 

1.0973 

15.000 

0.5701 

0.4031 

75.000 

1.6651 

1.1774 

20.000 

0. 6680 

0.4723 

80.000 

1. 7941 

1.2686 

25. 000 

0. 7585 

0.5363 

85. 000 

1.9479 

1.  3774 

30.000 

0. 844G 

0. 5972 

8G.46G 

2. 0000 

1.4142 

35.000 

0.9282 

0.65G3 

90.000 

2.1460 

1.5175 

39.3*17 

1.0000 

0.7071 

95.000 

2.4477 

1. 7308 

40.000 

1.0108 

0.7147 

95.450 

2.  4860 

1.7579 

45.  000 

1.0955 

0.7732 

98. 000 

2.  7971 

1.9778 

50.000 

1.1774 

0.8325 

98. 1G8 

2.8284 

2.0000 

54.40G 

1.2533 

0.  88G2 

98.889 

3. 0000 

2.1213 

55. 000 

1.2G37 

0.893G 

99.000 

3.0348 

2.1460 

GO. 000 

1.  3537 

0.9572 

99. 730 

3.4393 

2.4320 

The  probability  ellipse  that  contains  p-percent  of  the  wind  vectors 
expressed  in  the  most  general  form  is  the  conic  defined  by 

AX2  +  BXY  +  CY2  +  DX  +  EY  +  F  =  0  ,  (1 

where 

A  =  Oy2 
B  =  *2poxoy 

C*°x2 

D  =  2o  o  pY  -  2o  2X  =  -  (B?  +  2AX) 

x  y  y 

E  =  2o  o  pX  -  2a2Y  =  -  (BX  +  2C?) 
x  y  x 

F  =  AX2  +  CY2+  BXY  -  AC  (1  -  p2) 


and 


For  graphical  presentations,  the  range  of  the  variable  is  important  in 
order  to  arrange  the  scale.  The  largest  and  smallest  values  of  X  and  Y  for 
a  given  probability  ellipse,  p,  are  given  by 


where,  as,  before,  Xfi 


Although  there  are  several  approaches  to  graphing  the  probability  ellipses, 
the  following  procedure  is  advantageous  for  electronic  computer  plotting.  In 
establishing  the  computer  plotting  program,  the  sample  estimates  for  X,Y,  , 

»  ,  and  p  are  constants  in  equation  (13).  The  user  makes  the  choice  of  prob¬ 
ability  ellipses  desired.  Thus,  p  in  equation  (12)  is  programmed  as  a  parameter. 
The  largest  and  smallest  values  for  X  and  Y  are  computed  by  equations  (14)  and 
(15)  for  the  largest  probability  ellipse  selected.  This  sets  the  graphical 
scale.  Values  of  X  within  the  range  of  "X  smallest"  to  "X  largest"  are  obtained 
by  incrementing  X  between  these  limits.  Using  the  quadratic  equation,  a  solution 
for_Y  of  equation  (13)  is  made  and  plotted  for  each  value  of  X.  The  centroid 
(X,Y)  for  the  family  of  probability  ellipses  is  plotted  as  a  point.  Labeling 
and  other  identification  complete  the  plotting  program. 

For  a  given  probability,  equation  (13)  defines  an  ellipse  that  contains 
p- percent  of  the  points  X,Y.  Since  the  entire  area  under  the  bivariate  normal 
density  function  [equation  (9)]  is  unity,  upon  integration  for  a  given  prob¬ 
ability  ellipse,  that  given  ellipse  contains  p-percent  of  the  total  area.  In 
the  wind  statistics,  p-percent  of  the  wind  vectors  fall  within  the  specified 
probability  ellipse.  From  this  point  of  view,  a  specified  probability  ellipse 
gives  the  joint  probability  that  p-percent  of  the  U-V  components  lie  within 
the  given  ellipse. 

When  a  2  =  a  2  =  a2  and  p  =  0  in  the  bivariate  normal  distribution,  the 
x  y 

probability  ellipses  of  equation  (13)  reduce  to  circles  whose  centers  are  at 
the  means  Y,Y.  The  radii  of  the  probability  circles  are  Oy^,  where 


Values  for  Xc  for  selected  probabilities,  p,  are  given  in  table  C. 

Because  this  function  is  simple,  it  can  easily  be  graphed  manually. 
However,  the  generalized  plotting  technique  for  electronic  computer  plotters, 
as  represented  by  equation  (13),  can  be  advantageously  used. 


C.3.  Derived  Distributions  for  Wind  Statistics 

In  this  subsection  the  probability  distribution  functions  and  sets  of 
equations  are  presented  to  derive  certain  probability  distribution  functions 
for  wind  statistics.  These  derived  probability  distributions  are: 

1)  The  conditional  distribution  of  wind  components 

2)  The  generalized  Rayleigh  distribution  for  windspeed 

3)  The  distribution  for  wind  direction 

4)  The  conditional  distribution  of  windspeed  given  a  wind  direction 
(wind  rose). 

The  required  five  statistical  parameters  for  these  derived  distributions 
for  wind  statistics  are  given  in  table  I. 

C.3.1  The  Conditional  Distribution  of  Wind  Components 

Given  that  two  random  variables  X  and  Y  are  bivariate  normally  distributed, 
the  conditional  distribution  f ( Y | X )  is  read  as  f ( Y )  given  X,  and  likewise  f(X|Y) 
is  read  as  f ( X )  given  Y.  The  conditional  probability  distribution  function 
F( Y J X )  has  the  mean  E( Y } X )  and  variance  a2,  .  where 


E(Y  |X*)  =  Y  + 


<i) 


The  conditional  standard  deviation  is 


°(y |x*)  =  °y  /l  -  P2  •  <20) 


By  interchanging  the  variables  and  parameters,  the  conditional  distribu¬ 
tion  function  for  F(X|Y*)  has  the  conditional  mean 
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E(X  )  Y*)  =  X  +  p^j  (Y*  -  Y) 


(21) 


conditional  variance 


°2(x|y*)  =  ox2  (1  -  P2)  •  (22) 

and  conditional  standard  deviation 

°Cx|y*)  ■  °x  ^  -  »2  ■  (23) 


The  preceding  conditional  probability  distribution  functions  are  univariate 
normal  distributions  for  a  (fixed)  given  value  for  one  of  the  bivariate  normal 
variables.  Thus,  the  t-values  given  in  table  B  are  applicable  for  conditional 
probability  statements.  For  example. 


F(Y  |X*)  =  E(Y  | X *)  +  to(y|x*} 


(24) 


For  t  =  1.6449  there  is  a  95  percent  chance  that  Y  is  less  than  or  equal  to 
Y  +  1.6449  °(y|x*)  given  that  X  =  X*.  In  symbols  this  statement  reads 


p{  Y  <E(Y|X*)  +  1.6449  o(yjx#)  |X  =  X*J  =  0.9500  .  (25) 

Interval  probability  statements  can  also  be  made;  namely, 

Pj  Yx  =  E(Y|X*)  -  to(y|x.}  <  Y  <  Y2  =  E(Y|X*)  +  toy  |X  =  X*  J 

where  X*  can  take  on  any  fixed  value  of  X,  but  a  convenient  arrangement  is  to 
let  X*  =  X  ±  tox. 

The  close  connection  of  the  regression  function  of  Y  on  X  to  the  conditional 
mean  for  the  bivariate  normal  distribution  is  noted;  namely. 


(26) 


y  =  v*pg 

Similarly,  the  regression  function  of  X  on  Y  is 

X  =  X  +  p  ^  ^  (Y  -  Y)  . 

These  are  linear  functions  and  express  the  same  results  as  would  be  obtained 
from  a  least-squares  regression  line. 

C.3.2.  The  Generalized  Rayleigh  Distribution  for  Windspeed 

If  two  random  variables,  X  and  Y,  are  bivariate  normally  distributed,  then 
the  probability  distribution  for  the  modulus,  R,  can  be  derived  in  terms  of  the 
five  parameters  that  define  the  bivariate  normal  distribution. 


(28) 


The  distribution  of  R  so  derived  is  called  a  generalized  Rayleigh  distribution 
because  there  are  no  restrictions  on  the  parameters.  For  applications  to  the 
RRA,  the  variable  R  is  recognized  as  windspeed  or  the  modulus  of  the  wind  vector. 

The  probability  density  function  for  R  is  expressed  as 


f(R)  =  aQRe  -alR< 


Va2R  >  10(a3R) 


QQ 

+  Ij^agR^)  *2k^a3R^  008 


R  >  0 


(29) 


The  functions,  Iq(')»  I^'),  and  I * )  are  the  modified  Bessel  functions  of 
the  first  kind  for  zero  order,  kth  order,  and  2kth  order.  The  coefficients 
are 
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where  o&2  and  a^2  are  the  rotated  variances  to  produce  zero  correlation  between 
X  and  Y.  c?a  and  crb  are  the  positive  and  negative  roots'*  of  the  expression 


o 


2 


_  1 

(+,-)  2 


„  2/  .2  2.. 

°y  >  •  4ox  °y  41  - 


> 


ai  =  (0X2  +  Oy2)'4*1  -  P2)  ax20y2 


a 


2 


2/  .  2 

-  oy  )  +  ip  ox 

4(1  -  p2)  ax2oy 


2  2l  1/2 
°y  J 


J!  This  computational  form  is  obtained  from  the  determinant 


°X  '  K  °x°yp 


o  o  o 
x  yK 


°y2  "  K 


where  K  is  o2^+  and  and  are  analogous  to  the  standard  deviation  of 
the  major  and  minor  axes  of  the  bivariate  normal  probability  ellipse . 
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and 


tan  ip 


Since  this  density  function  cannot  be  integrated  in  closed  form  from  zero  to  R, 
numerical  integration  is  used  to  obtain  practical  results  for  the  probability 
distribution  function;  i.e.. 


R  + 

F(R)  =  /  f(R)  dR 

0  (30) 


A  number  of  special  cases  can  be  obtained  from  the  general  Rayleigh  distribution 
[equation  (29)],  the  simplest  of  which  is  to  let  a  =  a  =  a  and  X  =  Y  =  0  with 
independent  variables  X  and  Y.  This  gives  x  y 


f(R)  =  e~R2/2a2 

a1 


(31) 


which  is  recognized  as  the  classical  Rayleigh  probability  density  function.  The 
density  function,  equation  (31),  can  be  integrated  in  closed  form  over  any  range 
of  the  variable  R.  Hence,  the  probability  distribution  function,  F(R),  for 
equation  (31)  is 


F(R)  =  1  -exp 


(32) 


C.3.3.  The  Derived  Distribution  *>f  Wind  Direction 

Considering  the  wind  as  a  vector  quantity  and  bivariate  normally  distri¬ 
buted,  the  wind  direction  can  be  derived.  This  is  done  by  first  writing  the 
bivariate  normal  probability  density  function  in  polar  coordinates  whose 
variables  are 
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r  s  /  +  y  is  the  modulus  of  the  vector  or  speed,  and  0  is  the  direction  of 
the  vector.  After  integrating  g(r,e)  over  r  =  0  to  ®,  the  probability  density 
function  of  0  is 


2.  This  expression ,  equation  (33) ,  in  Smith  (1976)  is  given  with  respect  to 
the  mathematical  convention  for  a  vector  direction. 


O  o 

where  a  ,  b,  c  ,  and  are  as  previously  defined  in  equation  (33)  and 


»(&!-*<»•  7=  / 


2  dt 


is  taken  from  tables  of  normal  distribution  functions  or  made  available  through 
a  computer  subroutine. 

If  desired,  equation  (34)  can  be  integrated  numerically  over  a  chosen  range 
of  o  to  obtain  the  probability  that  the  vector  direction  will  lie  within  the 
chosen  range;  i.e.. 


F(0)  = 


g(0)  d0 


(35) 


One  application  may  be  to  obtain  the  probability  that  the  wind  will  flow  from 
a  given  quadrant  or  sector  as,  for  example,  onshore. 

C.3.4.  The  Derived  Conditional  Distribution  of  Windspeed  Given  the  Wind 
Direction  (Wind  Rose) 

Continuing  with  the  considerations  in  section  C.3.3.  of  this  chapter,  the 
conditional  probability  density  function  (pdf)  for  windspeed,  r,  given  a  speci¬ 
fied  value  for  the  wind  direction,  0,  can  be  expressed  as 


f(r  |  0) 


(36) 


where  the  coefficients,  a^  and  b  and  the  function  $  are  as  previously 
defined  in  equation  (33)  and  in  equation  (34). 

From  equation  (36)  the  mode  (most  frequent  value)  of  the  conditional 
windspeed  given  a  specified  value  of  the  wind  direction  is  the  positive 
solution  of  the  quadratic  equation, 


a2  r2  -  br  -  1=0 


(37) 
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The  locus  of  the  conditional  modal  values  of  windspeed  when  plotted  in  polar 
form  versus  the  given  wind  directions  forms  an  ellipse. 

The  noncentral  moment  for  equation  (36)  is  expressed  as 

00 

,;  =/  rn  f(r  1 6)  dr  .  (39) 

0 


Now  the  first  noncentral  moment  is  identical  to  the  first  central  moment  or 
the  expected  value,  E  (r|e).  The  integration  of  equation  (39)  for  the  first 
moment  is  sufficiently  simple  to  yield  practical  computations  and  can  be 
expressed  as 


Hence,  equation  (40)  gives  the  conditional  mean  value  of  the  windspeed  given 
a  specified  value  for  the  wind  direction. 

The  integration  of  equation  (36)  for  the  limits  r  =  0  to  r  =  r*  gives  the 
probability  t^at  the  conditional  windspeed  is  <_  r*  given  a  value  for  the  wind 
direction,  e.  This  conditional  probability  distribution  (PDF)  can  be  written 
as 


where 

r.  *  [a  r* '  (»)] 
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By  definition,  equation  (41)  is  an  expression  for  a  "wind  rose."  Empirical 
wind  rose  statistics  are  often  tabulated  or  graphically  illustrated  giving 
the  frequency  that  the  windspeed  is  not  exceeded  for  those  windspeed  values 
that  lie  within  assigned  class  intervals  of  the  wind  direction.  After  eval¬ 
uation  of  equation  (41)  for  various  values  of  windspeed,  r*,  and  the  given 
wind  directions,  e,  interpolations  can  be  performed  to  obtain  various  per¬ 
centile  values  of  the  conditional  windspeed. 

For  the  special  case  when  b  In  equation  (33)  equals  zero  (i.e.,  for 

x  3  y  -  0),  the  conditional  modal  values  of  windspeeds  [equation  (38)],  the 
conditional  mean  values  of  windspeeds  [equation  (40)],  and  the  fixed  condi¬ 
tional  percentile  values  of  windspeeds  [interpolated  from  evaluations  of 
equation  (41)],  when  plotted  in  polar  form  versus  the  given  wind  directions, 
produce  a  family  of  ellipses. 

For  the  special  case  when  7  =  y  =  0,  equation  (36)  reduces  to  the  fol¬ 
lowing  simple  case: 


Pr 


lr- 


6=0 


}■ 


(42) 


There  is  a  special  significance  of  equation  (42)  when  related  to  the  bivar¬ 
iate  normal  probability  distribution.  If  r*  and  e  are  measured  from  the 
centroid  of  the  probability  ellipse,  then  the  probability  that  r  _<  r*  is  the 
same  as  the  given  probability  ellipse.  Further,  solving  equation  (42)  for 
r*,  gives 


a  /  - 2  In  (1  -  P) 


(43) 


If  a  probability  ellipse  P  is  chosen,  equation  (42)  gives  the  distance  of 
r  along  any  e  from  the  centroid  of  the  ellipse  to  the  intercept  of  the 
specified  probability  ellipse.  If  there  is  an  interest  in  conditional  prob¬ 
ability  of  winds  for  a  given  9  relative  to  the  monthly  means,  equation  (43) 
is  applicable.  If  it  is  desired  to  find  the  magnitude  of  the  wind  along  any 
e  relative  to  the  monthly  mean  to  the  intercept  of  a  given  probability  ellipse, 
equation  (43)  is  applicable. 

D.  Statistical  Parameters  With  Respect  To  Any  Orthogonal  Axes 

The  five  wind  statistical  parameters  presented  in  table  I  are  given  with 
respect  to  the  standard  meteorological  coordinate  system;  i.e.,  these  param¬ 
eters  are  for  the  U  and  V  components.  For  many  aerospace  vehicles  and  range 
applications,  there  is  a  need  for  wind  statistics  with  respect  to  orthogonal 
axes  other  than  west  to  east  and  south  to  north.  For  example,  it  may  be 
required  to  present  wind  statistics  with  respect  to  a  flight  azimuth  of  an 
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aerospace  vehicle  whose  flight  azimuth  is  a  degrees  from  true  north  measured 
in  a  clockwise  direction.  The  following  sets  of  equations  are  presented  to 
compute  the  five  parameters  for  the  new  coordinate  axes  rotated  a  degrees 
clockwise  from  true  north. 

a.  Rotation  of  the  means  through  n  degrees: 

XQ  =  X  cos  (90  -  a)  +  X  sin  (90  -  a)  (44) 


Ya  =  Y  cos  (90  -  a)  -  X  sin  (90  -  a) 


b.  Rotation  of  the  variances  through  a  degrees: 


o  2  =  o  2  cos2  (90  -  a)  +  a  2  sin2  (90  -  a) 
a  y 


+  2pox°y  cos  (90  "  a>  sin  (9°  ~  a> 


2  2  2  „ _  2  2 

av  =  Gv  cos  (9°  “  a)  +  ovz  shi  (90  -  a) 


-  2P<Jxoy  cos  (90  -  a)  sin  (90  -  a)  .  (47) 


c.  Rotation  of  the  linear  correlation  coefficient  through  a  degrees: 


cov  (X ,Y) 

p  =  — - - i 

a  oo 
x  y 
a  Ja 


where  cov  (X,Y)a  is  the  rotated  covariance. 


cov  (X,Y)a  =  cov  (X , Y)  [cos2  (90  -  a)  -  sin2  (90  -  a)} 


+  cos  (90  -  a)  sin  (90  -  a)  (o  2  -  o  2) 

y  * 
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cov  (X ,Y)  =  P0x0y 


By  using  these  rotational  equations,  the  bivariate  normal  distribution  with 
respect  to  any  desired  rotated  coordinates  can  be  obtained  from  sample  esti¬ 
mates  that  have  been  computed  with  respect  to  a  specific  axis.  The  marginal 
distributions  after  rotation  are  also  normally  (univariate)  distributed. 
Using  the  rotational  equations  greatly  reduces  computational  efforts  for 
applications  requiring  statistics  with  respect  to  several  coordinate  axes. 

Appendix  A  presents  some  illustrative  examples  for  the  wind  statistics 
of  the  specific  RRA. 


CHAPTER  III.  STATISTICS  OF  THERMODYNAMICS  QUANTITIES 

AND  MODELS 

A.  General  Considerations 
A.l.  Objectives 

The  objective  inherent  in  developing  the  thermodynamic  section  of  the  RRA 
was  to  describe  the  thermodynamic  characteristics  of  the  atmosphere  using  a 
minimum  of  data  tabulations.  A  set  of  parameters  was  selected  which,  together, 
thermodynamically  describe  the  climatological  state  of  the  atmosphere.  These 
parameters  are  the  pressure,  temperature,  density,  dewpoint,  virtual  temper¬ 
ature,  and  water  vapor  pressure.  Used  together,  these  parameters  permit  the 
calculation  of  a  large  number  of  derived  quantities.  (Symbols  used  in  the 
calculations  in  this  chapter  are  summarized  in  table  D.)  Some  of  these 
quantities,  such  as  the  speed  of  sound,  are  dealt  with  in  section  III.E. 

The  probability  distribution  of  each  of  the  six  thermodynamic  RRA  param¬ 
eters  is  described  by  its  mean  value,  its  standard  deviation,  and  its  skewness. 
Several  of  these  parameters  (temperature,  pressure,  dewpoint  and  density)  have 
probability  distributions  that  are  close  to  a  univariate  normal  distribution; 
the  others  do  not.  The  skewness  parameter  gives  an  estimate  of  the  asymmetrical 
departures  of  a  probability  distribution. 

Hydrostatically  modeled  mean  values  of  pressure  and  density  were  calculated 
(table  IV),  so  that  users  may  determine  the  departure  of  the  actual  climatolo¬ 
gical  values  of  these  parameters  from  hydrostatic  conditions.  This  was  done  by 
hydrostatically  integrating  the  pressure  from  the  lowest  RRA  data  level  to  the 
termination  altitude  of  the  particular  RRA. 

A. 2.  Data  Quality  Control 

Data  limits  derived  from  the  following  parameters  were  used  to  screen  the 
thermodynamic  portion  of  the  RRA  data  base:  temperature,  pressure,  dewpoint 
(for  the  0-  to  30-km  portion  only),  and  density  (for  the  30-  to  70-km  portion 
only).  These  limits  were  set  to  plus  and  minus  six  standard  deviations  from 
the  mean  values  of  each  of  these  quantities.  These  limits  were  used  to  screen 
the  thermodynamic  portion  of  the  RRA  data  base,  according  to  the  procedures 
described  in  section  I.C.  The  data  base  used  to  generate  the  thermodynamic 
portion  of  the  RRA  (tables  I,  II,  and  IV)  was  considered  to  be  free  from 
errors  under  the  following  conditions: 

a)  The  skewness  values  of  the  pressure  and  temperature  were  between 
-2.5  and  2.5  at  all  data  levels. 

b)  The  skewness  values  of  the  density  were  between  -3.5  and  3.5  at 
data  levels  between  0  and  30  km. 

c)  The  skewness  values  of  the  density  were  between  -3.0  and  3.0  at 
data  levels  between  30  and  70  km. 

d)  The  skewness  values  of  the  dewpoint  were  between  -2.5  and  2.5  at 
all  data  levels  with  more  than  10  data  values. 


TABLE  0.  LIST  OF  SYMBOLS  USED  IN  CHAPTER  III 


Speed  of  sound 
Collision  diameter 
Vapor  pressure 
Gravity  at  latitude  <p 
Geopotential  height 

Geopotential  height  at  a  mandatory  radiosonde  data  level 

Geopotential  height  at  a  significant  radiosonde  data  level 

Coefficient  of  thermal  conductivity 

Mean  free  path  length 

Mean  molecular  weight  of  air  at  sea  level 

Annual  or  monthly  third  moment  of  quantity  Q 

Refractive  modulus 

Refractive  index 

Avogadro's  constant 

Number  of  values  of  quantity  Q 

Pressure 

Pressure  at  a  mandatory  radiosonde  data  level 

Pressure  at  a  significant  radiosonde  data  level 

Hydrostatically  integrated  mean  monthly  or  annual  pressure 

Any  tabulated  RRA  quantity 

Universal  gas  constant 

Specific  gas  constant  of  dry  air 

Parameters  used  in  converting  z  to  h  and  vice  versa 


TABLE  0.  (concluded) 


-  Sutherland's  constant ,  used  in  the  calculation  of  dynamic 
viscosity 

Temperature 

Dew  point 

Virtual  temperature 

Virtual  temperature  at  a  mandatory  radiosonde  data  level 
Virtual  temperature  at  a  significant  radiosonde  data  level 
Mean  air  particle  speed 
Mean  collision  frequency 

Parameter  used  in  the  hydrostatic  interpolation  of  pressure 
and  density 

Geometric  altitude 

Wavelength 

Skewness  of  quantity  Q 

Constant  used  in  the  equation  for  viscosity 

Ratio  of  specific  heat  at  constant  pressure  to  specific  heat  at 
constant  volume 

Kinematic  coefficient  of  viscosity 
Dynamic  coefficient  of  viscosity 
Density 

Mean  monthly  or  annual  density  derived  from  pressure  height 
Standard  deviation  of  the  quantity  Q 


A.  3  Limitation  of  Thermodynamic  Statistics 

The  correlation  coefficients  between  the  thermodynamic  quantities  and 
the  moisture-related  quantities  were  not  calculated  at  discrete  altitudes, 
nor  were  any  of  the  correlations  between  altitudes.  Therefore,  valid  sta¬ 
tistical  dispersion  models  that  require  the  relationship  between  two  or  more 
of  these  quantities  at  the  same  altitude  or  between  altitudes  cannot  be 
derived.  Approximations  for  the  correlation  coefficients  between  pressure, 
virtual  temperature,  and  density  at  discrete  altitudes  may  be  obtained  from 
the  coefficients  of  variation  as  developed  by  Buell  (1970).  The  coefficient 
of  variation  is  the  standard  deviation  divided  by  the  mean.  The  mean  values 
and  the  standard  deviations  are  taken  from  table  II.  A  model  for  the  profile 
of  monthly  and  annual  mean  pressure,  virtual  temperature,  and  density  that  is 
in  excellent  agreement  with  the  respective  statistical  mean  values  is  given 
by  table  IV.  This  agreement  results  because  the  physical  relationships,  given 
by  the  hydrostatic  equation  and  the  equation  of  state,  were  used  to  derive 
table  IV.  When  only  the  monthly  or  annual  mean  values  for  pressure,  virtual 
temperature,  and  density  are  required,  it  is  recommended  that  table  IV  be  used. 

B.  Establishing  Data  Samples  at  the  Required  Altitude  Levels 

This  section  describes  the  computational  procedures  used  to  establish  data 
samples  of  the  thermodynamic  RRA  parameters  at  the  RRA  data  levels.  References 
are  cited  only  when  an  equation  given  is  one  of  many  available  in  the  literature 
or  when  an  equation  is  stated  in  an  unusual  form. 

B.l.  Conversion  of  Data  Recorded  in  Geopotential  Heights  to  Geometric  Altitude 

The  upper-air  rocketsonde  observations  used  to  obtain  the  table  values 
above  30  km  were  recorded  in  terms  of  geometric  altitude  and  can  be  interpo¬ 
lated  directly  to  the  altitude  intervals  shown  in  the  tables.  However,  the 
radiosonde  observations  used  to  obtain  the  tabular  values  below  30  km  were 
recorded  in  terms  of  geopotential  heights.  The  change  of  coordinates  from 
geopotential  heights  to  geometric  altitudes  (h  to  z)  is  accomplished  by  cal¬ 
culating  a  table  of  geopotential  heights  that  correspond  exactly  to  the 
geometric  altitudes  at  which  the  atmospheric  parameters  are  tabulated.  The 
radiosonde  observations  are  then  interpolated  to  these  geopotential  heights. 

The  relationship  used  to  calculate  geometric  altitude  from  geopotential 
height  is 

H  =  (r'z)/(r*z)  ,  (49) 


where 

r  =  gr*/9. 80665 


and 

r*  =  -2g^J(d g^I d zq) 
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is  the  sea-level  gravity  at  the  latitude  <j>  corresponding  to  the  proper  loca¬ 
tion.  This  value  is  given  by  (List,  1968) 


=  9.780356  (1  +  5.2885  x  io'3  sin2<}>  -  5.  9  x  io"6  sin2  (2<J>)).  (50) 


■= — —is  the  rate  of  change  of  gravity  at  the  sea  level.  This  quantity  is  given 

o  Z_ 


by  the  equation 

---l  -  -3. 085462  x  10'6  +  2.27  x  10'9  cos  (2<j>)  -  2  x  10’12  cos  (4<J>)  .  (51) 
3zo 

The  units  used  for  gravity  are  meters  per  square  second,  while  the  units  for 

3  S<j, 

- — are  per  square  second. 

c)  z 


The  resulting  table  of  values  of  H  obtained  by  using  even  increments  of 
2  in  equation  (49)  is  shown  in  table  IV  of  the  RRA.  The  values  of  H  above 
30  km  are  not  used  in  the  interpolation  of  original  data,  but  are  included 
for  the  convenience  of  the  user. 

B.2.  Calculations  on  the  Original  Rawinsonde  Data  Records 

It  was  necessary  to  interpolate  the  information  from  the  original  rawin¬ 
sonde  data  records  to  the  geometric  altitudes  specified  as  the  RRA  data  levels 
The  parameters  for  which  this  interpolation  was  required  were  the  temperature, 
dewpoint,  and  pressure.  The  other  parameters  were  calculated  from  the  inter¬ 
polated  values  at  each  RRA  data  level.  These  "derived"  parameters  were  the 
water  vapor  pressure,  density,  and  virtual  temperature. 

B.2.1.  Calculation  of  the  Geopotential  Height  at  Significant  Levels 

Two  somewhat  different  interpolation  procedures  were  used  to  obtain  data 
from  radiosonde  and  rocketsonde  observations  at  the  levels  shown  in  the  tables 
The  procedure  used  to  interpolate  radiosonde  observations  began  with  the  cal¬ 
culation  of  virtual  temperature  at  each  data  level  in  a  sounding.  The  virtual 
temperature  was  computed  by 


Tv  =  T/(l.  -  0.379  (e/p)) 


(52) 


where  Ty  and  T  are  in  degrees  Kelvin  and  e  and  p  are  in  millibars. 
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The  radiosonde  soundings  contain  a  mix  of  data  taken  at  "mandatory"  and 
"significant"  levels.  Pressure,  temperature,  and  dewpoint  information  was 
given  in  these  soundings  at  both  types  of  levels.  However,  geopotential 
height  information  was  only  given  at  the  mandatory  levels.  The  heights  at 
the  significant  levels  were  "filled  in"  (calculated)  hydrostatically  using 
pressure  and  temperature  data  from  these  levels.  This  procedure  permitted 
the  use  of  most  of  the  significant  level  data  in  the  calculation  of  the  RRA 
tables.  The  equation  used  for  this  process  was 

T  T 

H  =  H_  +  29.2712617  *  --  — — 5 1  ln(P  /P  )  ,  (53) 

s  m  z  s  m 


where  the  subscripts  s  and  m  denote  quantities  at  significant  and  mandatory 
levels.  This  equation  was  not  used  if  the  difference  between  two  adjacent 
mandatory  levels  was  greater  than  200  mb.  All  soundings  with  such  data  gaps 
were  rejected  for  use  in  compiling  the  RRA. 

B.2.2.  Temperature 

Radiosonde  temperatures  were  interpolated  logarithmically  with  respect  to 
pressure  using  the  equation 


T 


+  <V  V 


lnp  -  lnpL 

‘"Pu  ■  lnpL 


(54) 


where  the  subscripts  U  and  L  indicate  values  at  the  nearest  data  levels 
in  the  actual  sounding  above  and  below  the  interpolated  level. 

B.2.3.  Pressure 

The  pressure  values  in  each  radiosonde  sounding  were  interpolated  to 
the  RRA  data  leva’s  using  the  equation 

_/  HL  -  Hu  ^ 

p  PL  exPl  29.2712617  (0.5)  (T  +  T  )  /  (55) 

VU  VL 


where  the  subscript  L  indicates  virtual  temperature,  geopotential  height, 
and  pressure  values  at  the  data  level  below  and  closest  to  the  level  at 
which  data  were  required. 

B.2.4.  Dewpoint  Temperature 

Dewpoint  values  were  interpolated  logarithmically  with  respect  to  pressure 
using  the  equation 

/lnp  -  lnp  \ 

Td  -  TdU  +  (TdL"  TdU}  (lnp”-  lnp "  j  *  (56) 
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The  subscripts  U  and  L  indicate  data  at  the  nearest  upper  and  lower  data  levels 
in  a  sounding. 

B.2.5.  Derived  Water  Vapor  Pressure 

The  water  vapor  pressure  was  calculated  from  the  interpolated  dewpoint 
values  at  the  RRA  data  levels  using  Teten's  approximation: 


7.5(Td  -  273.15)/(Td  -  35.86)  . 

e  =  6.11  mb  x  10  1  7) 

B.2.6.  Derived  Density 

The  density  values  derived  from  radiosonde  observations  were  calculated 
at  the  RRA  data  levels  using  the  equation 


p  =  348.36787  p/Tv 


(58) 


B.2.7.  Derived  Virtual  Temperature 

The  virtual  temperature  values  were  calculated  at  the  RRA  data  levels  for 
each  sounding  using  the  equation 

Ty  =  17(1  -  0. 379(e/p) )  ,  (59) 

where  Ty  and  T  are  in  degrees  Kelvin,  and  p  and  e  are  the  pressure  and  vapor  pres 
sure,  respectively,  in  millibars. 

B.3.  Calculations  on  the  Original  Rocketsonde  Data  Records 

The  rocketsonde  data  records  used  to  calculate  the  RRA  table  values  above 
30  km  were  given  in  terms  of  geometric  altitude.  For  this  reason,  slightly 
different  calculations  were  required  to  convert  the  recorded  data  values  to 
values  at  the  RRA  data  levels.  The  pressure,  temperature,  and  density  were 
all  interpolated  to  the  RRA  data  levels;  moisture-related  parameters  (virtual 
temperature,  water  vapor  pressure,  and  dewpoint)  were  not  calculated,  since 
atmospheric  moisture  at  altitudes  above  30  km  was  considered  to  be  negligible. 

No  interpolation  was  done  across  gaps  in  the  pressure  or  temperature  data 
within  a  sounding  larger  than  7,000  m.  Data  values  at  the  RRA  levels  within 
such  a  gap  were  set  to  missing. 

B.3.1.  Temperature 

Rocketsonde  temperatures  were  interpolated  linearly  with  respect  to  geo¬ 
metric  altitude  using  the  equation 
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T  =  Ty  +  (tl  V 


(60) 


where  the  subscripts  U  and  L  indicate  values  at  the  nearest  data  level  in  the 
actual  sounding  above  and  below  the  interpolated  level. 

B.3.2.  Pressure 

The  pressure  values  in  each  rocketsonde  sounding  were  interpolated  to  the 
RRA  data  levels  using  the  equation 


where 


=  -  TyU  +  TvL 

1  v  2 


B.3.3.  Density 

Rocketsonde  density  values  were  interpolated  using  the  equation 


p  =  PL  exp 


(Z  -  ZL> 
Tf*~  fP" 


(62) 


where  W  is  specified  in  section  1 1 1. B.3.2. 

C.  Computation  of  Statistical  Parameters  for  Tables  II  and  III 

A  three-step  procedure  was  used  for  computing  the  monthly  and  annual 
means,  standard  deviations,  and  skewness  values  from  the  data  values  at  the 
RRA  data  levels.  Initially,  certain  statistical  sums  were  calculated  and 
stored  as  the  soundings  in  the  data  base  were  processed.  These  sums  were 
then  used  to  calculate  the  monthly  statistics  given  in  the  RRA  tables.  The 
annual  statistics  were  then  calculated  from  these  stored  sums  and  the  monthly 
statistics . 

C.l.  Stored  Statistical  Sums 
The  sums  calculated  were 


* 

* 
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£Q.  £Q2.  and  £q3 


where  Q  is  any  one  of  the  quantities  given  in  the  thermodynamic  part  of  the 
RRA. 

C.2.  Calculation  of  the  Monthly  Statistics 
C.2.1 .  Monthly  Means 

The  mean  monthly  values  of  the  thermodynamic  RRA  quantities  were  calcu¬ 
lated  using  the  equation 

Q  =  Eq,nQ  1 

where  Ng  is  the  number  of  observed  values  of  the  quantity  Q  for  a  given  month. 
C.2.2.  Monthly  Standard  Deviations 

The  monthly  standard  deviations  of  the  thermodynamic  RRA  quantities  were 
calculated  using  the  equation 


(NqE'Q2)  -  (IQ) 2 


Q 


(Nq  -  1) 


(63) 


C.2.3.  Monthly  Skewness  Values 

The  monthly  skewness  values  of  the  windspeed  and  of  the  thermodynamic  RRA 
quantities  were  calculated  using  the  equation 

M3 

sv  —  ^ 


where  M3n  is  the  third  moment  of  the  quantity  Q,  on  is  its  standard  deviation. 


C.3.  Calculation  of  the  Annual  Statistics 

Equations  (63)  and  (64),  used  to  calculate  the  monthly  values  of  the 
standard  deviations  and  skewness  values,  involve  taking  the  differences 
between  two  pairs  of  large  sums  containing  Q2  and  Q3,  where  Q  is  the  ther¬ 
modynamic  RRA  quantity.  Using  these  equations  to  compute  the  annual  statis¬ 
tics  would  have  resulted  in  a  substantial  loss  of  precision,  as  these  sums 
become  larger  by  several  orders  of  magnitude  in  such  a  case.  This  problem 
was  avoided  by  calculating  the  annual  means,  standard  deviations,  and  skew¬ 
ness  values  from  the  monthly  statistics. 

C.3.1  Annual  Mean  Values 

The  annual  mean  values  of  the  thermodynamic  RRA  quantities  were  calculated 
using  the  equation 

QANN  =  QA/NQ 


where  is  the  total  of  all  observed  values  of  Q  and  Nq  is  the  total  number  of 
observations  of  Q. 

C.3.2.  Annual  Standard  Deviations 

The  annual  standard  deviations  of  the  thermodynamic  RRA  quantities  were 
calculated  using  the  equation 


where  ^  =  the  number  of  data  values  for  Q  in  month  i  (i  =1  to  12),  Qi  =  the 
monthly  mean  of  Q,  and  Op..  =  the  standard  deviation  of  quantity  Q  in  month  i. 

C.3.3.  Annual  Skewness  Values 

The  annual  skewness  values  of  the  thermodynamic  RRA  quantities  were  calcu¬ 
lated  using  the  equation 


12 

12 

M3Qann  =  k  X 
1=1 

<NQ! 

M3Qi)+N«ANN  iS 

-  nji  1 

12 

v 

/K  n  .  ^ANN 

'  Qi  '  Mb 

^ANN 

*»ANN 

k 

12 


^ANN  12 


NQ 


ANN 


2  v  3 


(NQi  °Qi  >  +  2§ANN 


(66) 


where  =  the  third  moment  about  the  mean  of  quantity  Q  in  month  i  and 
M3QaNN  =  the  annual  third  moment  about  the  mean  of  the  quantity  Q. 

D.  Derived  Monthly  Mean  and  Annual  Mean  Model  Atmospheres 

A  set  of  modeled  monthly  mean  and  annual  mean  hydrostatic  values  of  pres¬ 
sure  and  density  was  calculated  from  the  lowest  RRA  data  level  (0  km,  mean  sea 
level)  upwards  to  30  km,  and  from  30  km  upwards  to  70  km.  The  integration  from 
0  to  30  km  was  computed  independently  of  the  integration  from  30  to  70  km  because 
of  the  difference  in  data  sources.  The  two  different  values  for  30  km  are  pro¬ 
vided  for  comparison.  When  30-km  data  are  required,  the  values  given  in  the 
0-  to  30-km  table  should  be  used.  These  hydrostatically  modeled  mean  values, 
which  are  given  in  table  IV,  are  useful  as  a  check  on  the  validity  of  the  pres¬ 
sure  and  density  values  given  in  table  II.  In  most  cases,  the  values  in  tables 
II  and  IV  for  any  given  data  level  are  within  1  percent  of  each  other.  The 
hydrostatic  pressure  values  in  table  IV  were  calculated  using  the  equation 

/  0.034162  (Hj  -  HQ) 

P1  =  p0  exp  (  '  0.  5  (T  +T  ) 

\  V1  V0 


where  -  Hq  is  in  meters  and  a  "0"  subscript  refers  to  values  at  the  RRA  data 

level  immediately  below  the  level  being  checked,  pg  at  the  lowest  data  level  is 

set  equal  to  the  RRA  mean  pressure;  p^ ,  calculated  for  the  next  highest  data 

level,  is  taken  as  pQ  for  the  level  above  that.  This  process  is  repeated  for 

all  the  other  RRA  data  levels.  The  hydrostatic  density  corresponding  to  the 
hydrostatic  pressures  is  calculated  from  these  pressures  and  the  RRA  virtual 
temperature  values  using  the  formula 


P„  =  348.36786  P„/T 


(68) 


where  ^  and  P^  are  the  hydrostatic  density  and  pressure  shown  in  table  IV  of 
the  RRA. 

E.  Thermodynamic  Quantities  Derivable  from  the  Basic  Tables 

Several  other  quantities  can  be  calculated  from  the  statistics  listed  in 
tables  I  and  II.  Primary  physical  constants  used  in  these  calculations  are 
listed  in  table  E.  The  equations  given  in  this  section  can  be  used  to  calcu¬ 
late  the  approximate  mean  values  of  these  quantities  at  each  RRA  data  level. 
It  is  not  possible  to  infer  or  derive  any  information  concerning  the  standard 
deviation  or  skewness  values  of  these  quantities  from  the  data  in  tables  II 
and  III  of  the  RRA. 

E.l.  Mean  Air  Particle  Speed 

The  mean  air  particle  speed,  V,  is  the  arithmetic  average  of  the  speeds 
of  all  air  particles  in  the  volume  element  being  considered.  For  a  valid 
average  to  occur,  there  must  be  a  sufficient  number  of  particles  involved 
to  represent  mean  conditions.  The  equation  for  V  for  dry  air  is 

V  =  A  . 

/  tr  M 


A  computational  form  tor  dry  air,  using  tabulated  values,  is 


V  =  /  7  .3094  x  102  x  T  (meters  per  second)  , 


(69) 


(70) 


where  T  is  the  temperature  in  degrees  Kelvin  from  table  II.  Equation  (69), 
when  corrected  fol"  moist  air,  becomes 


V  = 


R'  Ty 


(71) 


The  computational  form  for  moist  air  is 


V 


/ 7.3094 


Ty  (meters  per  second) 


(72) 


where  Ty  is  the  virtual  temperature  in  degrees  Kelvin  from  table  III. 


TABLE  E.  LIST  OF  PRIMARY  PHYSICAL  CONSTANTS 


=  standard  atmospheric  pressure  at  sea  level 
1.013250  x  io5  Newton/m2  =  2116.22  lb/ft2 

standard  atmospheric  density  at  sea  level 
1.2250  kg/m3  =  0.076474  lb/ft3 

standard  temperature  at  sea  level  =  288.15  K  =  15. 0°C  —  59.0°F 

standard  gravity  at  sea  level  at  latitude  45°32,33n 
9.  80665  m/s2 

Sutherland’s  constant  used  in  calculation  of  dynamic  viscosity 

110.4  K 

ice-point  temperature  at  PQ  =  273.15  K 

constant  used  in  calculation  of  dynamic  viscosity 
1.458  x  10'6  kg/s  m  K* 

7.3025  x  10' 7  ]b/s  ft  R* 

ratio  of  specific  heat  of  air  at  constant  pressure  to  specific 
heat  of  air  at  constant  volume 

1.4 

mean  effective  collision  diameter  of  air  molecules 

-in  -  q 

3.65  x  10  m  =  1.1975  x  10  ft 

Avogadro's  constant 

6.022169  x  1026/kg  mol  =  2.73179  x  1026/lb  mol 
gas  constant  =  8.31432  J/mol  K 

2 

gas  constant  for  dry  air  =  2.  8704  x  10  J/kg  K 
molecular  weight  of  dry  air  =  28.966  g7mol 


The  mean  free  path,  L,  is  the  mean  value  of  the  distance  traveled  by  each 
neutral  air  particle  in  a  selected  air  parcel,  between  successive  collisions 
with  other  particles  in  that  parcel.  A  meaningful  average  requires  that  the 
selected  parcel  be  large  enough  to  contain  a  substantial  number  of  particles. 
The  equation  for  L  is  given  by 


where  C.  is  the  effective  collision  diameter  of  the  mean  air  molecules.  The 

d  -10 
1976  standard  atmosphere  value  of  3.65  x  10  is  valid  for  the  range  of  alti¬ 
tudes  in  the  RRA. 

A  computational  form  for  moist  air,  using  tabulated  values,  is 

L  =  2.335  x  10  7  |  (meters)  ,  (74) 

where  T  is  the  temperature  in  degrees  Kelvin  from  table  II  and  P  is  the  pres¬ 
sure  in  millibars  from  table  II. 

A  form  of  (73)  to  correct  L  for  moist  air  is 

L-(Q)^  . 

'  '»«  C„2  (75) 

The  computational  form  for  moist  air  is 

-7  Tv 

L  =  2.3325  x  10  -p  (meters)  ,  (76) 

where  Ty  is  the  virtual  temperature  in  degrees  Kelvin  from  table  III  and  P  is 
the  pressure  in  millibars  from  table  II. 

E.3.  Mean  Collision  Frequency 

The  mean  collision  frequency,  Vc,  is  considered  to  be  the  average  speed 

of  air  particles  contained  in  an  air  parcel,  divided  by  the  mean  free  path  of 
the  particles  inside  that  parcel.  Computationally  this  is  equivalent  to 


V  =  f-  (sec  ) 
C  L 


(77) 


To  determine  Vc  for  dry  air,  use  V  and  L  from  equations  (70)  and  (74). 
To  determine  Vc  for  moist  air,  use  V  and  L  from  equations  (72)  and  (76). 

E.4.  Speed  of  Sound 


The  expression  for  the  speed  of  sound,  C$,  in  meters  per  second  in  dry 
air,  is 


C  = 

s 


(78) 


To  compute  Cg  for  dry  air  from  tabulated  values,  use 


=  / 


4.0185  »  10  x  t  (meters  per  second) 


(79) 


where  T  is  the  temperature  in  degrees  Kelvin  from  table  II.  One  form  for  the 
speed  of  sound  in  moist  air  is 


C 

s 


(80) 


where  Ty  is  the  virtual  temperature  from  table  III.  A  computational  form  for 
moist  air  is 


C 

s 


/ 


4.0185  x  10Z  T 


(meters  per  second) 


(31) 


E.5.  Dynamic  Coefficient  of  Viscosity 

The  coefficient  of  dynamic  viscosity,  u,  is  defined  as  a  coefficient  of 
internal  friction  developed  where  gas  regions  move  adjacent  to  each  other  at 
different  velocities.  The  following  expression  is  taken  from  the  U.S.  Stand¬ 
ard  Atmosphere  (1976): 


The  computational  form  is 


(1.458  x  10~6)  T3y2 
T  +  110.4 


(kilograms  per  second 
per  meter)  , 


(83) 


where  T  is  the  temperature  in  degrees  Kelvin  from  table  II. 

E.6.  Kinematic  Coefficient  of  Viscosity 

The  kinematic  coefficient  of  viscosity,  designated  as  n,  is  defined  to  be 
the  ratio  of  the  dynamic  coefficient  of  viscosity  of  a  gas  to  its  density,  or 


n  •=  y/p 


(84) 


The  computational  form  is 


3  (square  meters 

n  =  1.0  *  10  \\J  p  per  second) 


(85) 


where  p  is  the  dynamic  coefficient  of  viscosity  from  equation  (83)  and  p  is 
the  density  in  grams  per  cubic  meter  from  table  II. 

E.7.  Coefficient  of  Thermal  Conductivity 

The  empirical  expression  used  for  the  coefficient  of  thermal  conductivity, 
designated  as  Kt,  is  given  in  the  1976  Standard  Atmosphere  as 


2.65019  x  1  Q~ 3  .  T372 
T  +  245.4  x  io'(12yT) 


(watts  per  meter 
per  degree  Kelvin)  . 


(86) 


where  T  is  in  degrees  Kelvin. 

E.8.  Refractive  Modulus  and  Refractive  Index 

The  refractive  modulus  or  refractivity  (Selby  and  McClatchey,  1975;  Smith 
and  Weintraub,  1953)  is  defined  as  N,  where 


N  =  (n  -  l)  •  io6 


(87) 


and  n  is  the  refractive  index. 
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For  microwave  frequencies  below  approximately  30  GHz  (equivalent  to  wave¬ 
lengths  above  1  cm),  N,  the  refractive  modulus,  is  given  by  the  empirical 
equation 

N  =  77. 6  ^r-  +  3.73  x  105  (dimensionless),  (88) 


where  E  and  P  are  in  millibars  and  T  and  are  in  degrees  Kelvin. 


The  following  expression  is  valid  for  the  visible  and  infrared  wavelengths 
shorter  than  approximately  30  ym  (0.03  mm). 


N  =  77.6  ij  +  0.584  (dimensionless)  ,  (89) 


where  \  is  the  wavelength  in  microns  and  T  is  in  degrees  Kelvin. 

The  expression  for  N  for  the  wavelength  from  0.03  mm  to  1  cm  is  an  extremely 
complex  function  of  wavelength. 


CHAPTER  IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

This  document  satisfies  the  technical  objectives  established  for  the  RRAC 
by  the  RCC  MG.  Upper  air  statistics  and  models  for  wind  and  thermodynamic  quan¬ 
tities  for  the  specific  site  have  been  derived  in  a  consistent  and  uniform  man¬ 
ner,  which  will  be  used  in  publications  for  all  other  assigned  site  locations. 
These  RRAs  represent  an  improvement  over  the  previously  published  RRAs  because 
of  the  availability  of  more  extensive  upper  air  data  bases  and  the  adaptation  of 
more  advanced  statistical  techniques.  A  statistical  measure  of  central  tendency 
(mean  values)  and  a  measure  of  dispersion  (standard  deviation  with  respect  to 
the  mean  values)  for  monthly  and  annual  reference  periods  have  been  tabulated 
for  all  variables  in  a  consistent  manner  from  data  bases  that  have  been  edited 
and  quality-controlled  in  the  same  manner.  Further,  a  statistical  measure  for 
symmetry  (skewness  coefficient  that  involves  the  third  statistical  moment)  has 
been  tabulated  for  all  variables  except  the  U  and  V  wind  components.  Even  with 
these  improvements,  the  user  of  these  RRAs  must  recognize  certain  limitations 
of  the  statistical  tabulations: 

1)  The  wind  profile  structure  with  respect  to  altitude  cannot  be  modeled 
from  the  RRA  statistics  because  the  interlevel  and  crosslevel  correlations  were 
not  computed. 

2)  The  profile  structure  with  respect  to  altitude  for  any  of  the  ther¬ 
modynamic  variables  or  any  quantities  derivable  from  these  variables  cannot 
be  modeled  because  the  prerequisite  correlations  were  not  computed.  However, 
the  profiles  of  monthly  and  annual  means  for  pressure,  virtual  temperature, 
and  density  are  in  agreement  (table  IV)  with  the  hydrostatic  equation  and 
the  equation  of  state. 

The  preceding  limitations  are  cited  to  prevent  a  misuse  of  the  RRAs. 

More  extensive  statistical  tabulations  were  beyond  the  scope  of  this  commit¬ 
tee's  task.  As  greater  insight  is  gained  through  usage  of  these  RRAs,  many 
adaptations  of  the  statistical  tabulations  for  specific  engineering  and  sci¬ 
entific  applications  are  envisioned. 

Recommendations 

It  is  recommended  than  the  wind  and  thermodynamic  statistical  tabulations 
and  attendant  models  contained  in  the  RRAs  be  used  as  a  standard  reference 
source,  as  may  be  appropriate,  by  the  ranges  and  range  users.  It  is  further 
recommended  that  the  respective  Range  Staff  Meteorologist  or  responsible  agency 
staff  member  be  consulted  for  the  applicability  of  the  RRAs  for  specific  engi¬ 
neering  applications. 
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Research  Laboratories.  Available  through  the  National  Technical  Information 
Service,  Washington,  D.C.,  1975. 

Smith,  E.K.;  and  Weintraub,  S.:  "The  Constants  in  the  Equation  for  Atmospheric 
Refractive  Index  at  Radio  Frequencies,"  Proceedings  of  the  Institute  of  Radio 
Engineers,  41,  8,  August  1953,  pp.  1035-1037. 

Smith,  O.E.:  NASA  TM  X-73319,  Vector  Wind  and  Vector  Wind  Shear  Models  at 
0-27  km  Altitude  for  Cape  Kennedy,  Florida,  and  Vandenberg  AFB,  California. 
Prepared  under  sponsorship  of  the  National  Aeronautics  and  Space  Adminis¬ 
tration.  Available  through  the  National  Technical  Information  Service, 
Washington,  D.C.,  July  1976. 

U.S.  Standard  Atmosphere,  1976.  Prepared  under  the  sponsorship  of  the 
National  Aeronautics  and  Space  Administration,  United  States  Air  Force,  and 
United  States  Weather  Bureau.  Available  through  U.S.  Government  Printing 
Office,  Washington,  D.C.,  October  1976. 

PREVIOUS  RANGE  REFERENCE  ATMOSPHERES  PUBLISHED  BY  IRIG 


Atlantic  Missile  Range  Reference  Atmosphere  for  Cape  Kennedy,  Florida  (Part  I), 
Document  104-63,  April  1 6,  1963.  (AD451780) 

White  Sands  Missile  Range  Reference  Atmosphere  (Part  I),  Document  104-63, 

June  28,  1964.  (AD451781) 

Fort  Churchill  Missile  Range  Reference  Atmosphere  for  Fort  Churchill,  Manitoba, 
Canada  (Part  I),  Document  104-63,  August  7,  1964.  (AD634727) 

Pacific  Missile  Range  Reference  Atmosphere  for  Eniwetok,  Marshall  Islands 
(Part  I),  Document  104-63,  September  1,  196,4.  (AD479264) 

Fort  Greely  Missile  Range  Reference  Atmosphere  (Part  I),  Document  104-63, 
October  6,  1964.  (AD634726) 

Eg! in  Gulf  Test  Range  Atmosphere  for  Eg! in  AFB,  Florida  (Part  I),  Document 
104-63,  January  25,  1965.  (AD472601) 

Pacific  Missile  Range  Atmosphere  for  Point  Arguello,  California  (Part  I),  Docu- 
ment  104-63,  April  1965.  (AD472602) 


ss.<. 


v»» 


Wallops  Island  Test  Range  Reference  Atmosphere  (Part  I),  Document  104-63, 


here  (Part  I),  Document  104-63, 


September  1971.  (AD751581 

White  Sands  Missile  Range  Reference  Atmosphere  (Part  II),  Document  104-63 
September  1 971 .  ( AD782654 ) 

Wallops  Island  Test  Range  Reference  Atmosphere  (Part  II),  Document  104-63 
September  1 971 .  (ADA040280) 

Fort  Greely  Missile  Range  Reference  Atmosphere  (Part  II),  Document  104-63 


September  1971 .  (ADA040281 

Edwards  Air  Force  Base  Reference  Atmosphere  (Part  I),  Document  104-63, 


September  1972.  (AD782651 

Kwaialein  Missile  Range  Reference  Atmosphere  for  Kwa.ialein,  Marshall  Islands 


(Part  I),  Document  104-63,  October  1974.  (ADA002664 

Pacific  Missile  Test  Center  Reference  Atmosphere  for  Point  Arguello,  Cali- 


ADA040279 


0-70  Km  Altitude,  Document  360-82,  December  1982.  CADI 23424 


Cape  Canaveral,  Florida,  Range  Reference  Atmosphere,  0-70  Km  Altitude,  Document 


361-83,  February  1983.  (ADA1 25553) 

Vandenberg  Air  Force  Base,  California,  Range  Reference  Atmosphere,  0-70  Km 


Altitude,  Document  362-83,  April  1983. 

Dugway,  Utah,  Range  Reference  Atmosphere,  0-30  Km  Altitude,  Document  363-83, 
June  1983. 

Wallops  Island,  Virginia,  Range  Reference  Atmosphere,  0-70  Km  Altitude,  Docu¬ 


ment  364-83,  July  1983. 
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White  Sands  Missile  Range,  New  Mexico,  Range  Reference  Atmosphere,  0-70  Km 
Altitude,  Document  365-83,  August  1983. 

Edwards  AFB,  California,  Range  Reference  Atmosphere,  0-70  Km  Altitude,  Docu¬ 
ment  366-83,  August  1983. 

Eglin  AFB,  Florida,  Range  Reference  Atmosphere,  0-30  Km  Altitude,  Document 
367-83,  September  1983. 

In  addition  to  the  documents  above  and  the  present  RRA  for  Taquac,  Guam  Island 
the  revised  series  will  include  RRAs  for  the  following  locations: 


Point  Mugu,  California 
Barking  Sands,  Hawaii 
Ascension  Island,  South  Atlantic 


Numerical  values  in  this  document  are  given  in  the  International  System 
of  Units  (SI,  Systeme  International  d'Unites).  The  values  in  parentheses  are 
equivalent  U.S.  Customary  Units,  which  are  English  units  adapted  for  use  by  the 
United  States  of  America.  The  SI  and  U.S.  Customary  Units  provided  in  table  F 
are  those  normally  used  for  measuring  and  reporting  atmospheric  data. 

By  definition,  the  following  fundamental  conversion  factors  are  exact: 

Type  U.S.  Customary  Units  Metric 

Length  1  U.S.  yard  (yd)  0.9144  meter  (m) 

Mass  1  avoirdupois  pound  (lb)  453.59237  gram  (g) 

Time  1  second  (s)  1  second  (s) 

Temperature  1  degree  Rankine  (°R)  9/5  degree  Kelvin  (K) 

To  aid  in  the  conversion  of  units,  conversion  factors  based  on  the  above 
fundamental  conversion  factors  are  given  in  table  F. 


TABLE  F.  FACTORS  FOR  CONVERSION  UNITS 


jtur.  wo;: 


mb  n.  75006  mnillr  <0*0 

in.  i32  I  >  25.40*  min  I  If  <(V’l » 

inniMf<0()  1.53322  mb 

in.ll*  071)  33.4639  mb 
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TABLE  |.  |  Mire  STATISTICAL  PARAt'CUS. 


STATION  - 

912170 

TACL'AC 

Z 

'ifA-i  u 

5.0.  U 

R1U.VI 

rCM  v 

S.O.  V 

MCAN  US 

S.O.  MS 

SKCH  MS 

NODS 

KM 

M/5 

M/5 

M/5 

M/5 

M/S 

.III 

-'•.Vi 

1.95 

.1641 

-1.19 

1 .82 

4.93 

i.et 

-.12 

375. 

1.000 

-9. 18 

4.80 

.2627 

-1.44 

4.16 

10.53 

3.95 

.05 

420. 

2.  COO 

-7.40 

5.25 

.2293 

-.66 

3.69 

8.96 

4.02 

.27 

422. 

2.CC0 

-6.34 

5.60 

.2652 

-.10 

3.80 

8.40 

3.93 

.33 

422. 

4.0C0 

-9.89 

5.52 

.2211 

.63 

3.9. 

8.11 

3.89 

.47 

4??. 

o.oco 

-9.80 

3.29 

.0429 

.71 

4.00 

8.03 

3.77 

.46 

423 . 

6  000 

-5.6? 

5.29 

-.0017 

1 .05 

4.61 

8.19 

3.83 

.39 

422. 

7.000 

-4.69 

5. 14 

-.0044 

.92 

4.45 

7.53 

3.51 

.62 

421 . 

9.000 

-3.65 

5.27 

-.0630 

.57 

4.29 

6.84 

3.60 

.77 

419. 

9.000 

-2.40 

5.08 

-.1159 

.42 

4.45 

6.29 

3.45 

.85 

4t8. 

10.000 

-1.43 

5.14 

-.1119 

.79 

4.51 

6.12 

3.45 

.99 

416. 

11.030 

-.71 

5.47 

-.0485 

1.59 

5.12 

6.70 

3.75 

.88 

416. 

lr?.  000 

-.70 

5.44 

-.0569 

2.87 

5.40 

7.25 

3.84 

.98 

413. 

I?. 000 

-1.12 

5.71 

- .0431 

4.1! 

3 .  l«J 

8.1  1 

4.16 

.9? 

413. 

I1*.  000 

-2.26 

5.76 

- . I 035 

5.27 

6.20 

9.12 

4.61 

.60 

413. 

>.  *0 . 000 

-4.4J 

5.87 

-.2696 

6.09 

6.78 

10.41 

5.36 

.52 

400. 

16.000 

-7.10 

6.01 

-.4467 

5.38 

5.82 

1C.  68 

5.94 

.65 

392. 

17.000 

-8.11 

5.38 

-.4980 

4. 15 

4.89 

10.26 

5.52 

.54 

349. 

10.000 

-6.84 

4.36 

-.2212 

1.62 

3.13 

7.97 

3.83 

.25 

350. 

n  ooo 

-5.29 

4.71 

-.0988 

.29 

2.42 

6.56 

3.62 

.42 

342. 

<■•0.000 

-4.10 

5.63 

.0171 

.  12 

1 .99 

6.13 

3.84 

.67 

341. 

?! .000 

-3.93 

6.8? 

-.1232 

.16 

2.  15 

6. 78 

4.54 

.63 

327. 

51 . COO 

-4.40 

7,85 

-.1116 

.24 

2.07 

7.86 

4.84 

.66 

320. 

23.000 

-4.80 

8.85 

-.0667 

.17 

2.31 

9.10 

4.88 

.84 

316. 

f'.OOO 

-5.30 

9.53 

-.0234 

.14 

2.45 

9.85 

5.25 

.81 

314. 

£•‘..000 

-4.84 

10.35 

-.0263 

-.02 

2.61 

10.09 

5.»t 

.85 

308. 

?9  coo 

-4.34 

11. 05 

-.0014 

.00 

2.99 

10.34 

6.54 

.93 

257. 

27 .  COO 

-2.34 

12.40 

.0460 

.10 

3.19 

11.19 

7.03 

1.13 

276. 

?<;.:;o 

-1.93 

13.49 

.07-0 

.41 

3.61 

12.33 

6  61 

I.C6 

256. 

29 .  COO 

-.e5 

14. C3 

.  1 255 

.57 

3.64 

13.07 

5.29 

.65 

229. 

20 . CCO 

-.23 

15.14 

.  1 55  3 

.86 

3.77 

14.07 

6.73 

.81 

169. 

TABLE 

I.  2 

HINO  STATISTICAL  parameters. 

FEBRUARY 

STATION 

-  912170 

TAGUAC 

Z 

MEAN  U 

S.O.  u 

R1U.V1 

ft  AN  V 

S.D.  V 

MCAN  US 

S.O.  US 

SKEW  HI*. 

NOBS 

KM 

M/S 

M/S 

M/S 

M/S 

M/S 

M/S 

.111 

-4.40 

1.79 

.1371 

-1 .25 

1 .81 

4.97 

1 .62 

.04 

341. 

1  .000 

-9.41 

3.91 

.2362 

-1.42 

3.75 

10.45 

3.27 

-.17 

383. 

2.000 

-7.21 

4.15 

.2816 

-1.09 

3.75 

8.47 

3.55 

.23 

384. 

3.  COO 

-4.98 

4.16 

.1451 

-.65 

4.01 

6.89 

3.34 

.89 

384. 

4.000 

-4.09 

4.03 

.1359 

-.55 

3.90 

6.22 

3.13 

.78 

383. 

5.000 

-4. 10 

4.13 

.0262 

-.74 

4.63 

6.69 

3.31 

.56 

383. 

6.000 

-3.62 

4.78 

.0099 

-.61 

5.  10 

6.92 

3.80 

.76 

382. 

7.000 

-3.25 

4.78 

- .0925 

-.88 

5.66 

7.27 

3.66 

.65 

377. 

8.000 

-2.15 

4.99 

-.2446 

-1,33 

5.64 

7.07 

3.62 

.81 

372. 

9.000 

-.40 

5.22 

-.31 14 

-1 .41 

5.29 

6.57 

3.76 

1.20 

372. 

10.000 

.91 

5.58 

-.1968 

-.64 

4.99 

6.65 

3.58 

1.28 

369. 

11.000 

1.89 

6.15 

-.1176 

.83 

4.83 

7.20 

3.68 

.77 

369. 

12.000 

2.36 

6.29 

-.0869 

3.00 

5.39 

S.12 

4.12 

.69 

367. 

13.000 

2.34 

6.55 

-.1273 

4.59 

5.90 

9.07 

4.67 

.67 

361. 

14. COO 

1.50 

6.72 

-.2161 

6.05 

6.26 

9.84 

5.11 

.61 

360. 

15.000 

-.24 

7.00 

-.3586 

5.63 

5.91 

9.61 

4.81 

.44 

355. 

16.000 

-2.97 

6.99 

-.5212 

4.13 

5.96 

8.95 

5.36 

.89 

344. 

17.000 

-4.24 

6.04 

-.4701 

2.53 

4.85 

7.81 

4.83 

.90 

328. 

18.000 

-3.79 

4.48 

-.2077 

.52 

3.09 

5.86 

3.  14 

.84 

330. 

19.000 

-2.79 

4.45 

-.0690 

-.18 

2.76 

5.09 

3.01 

1.00 

325. 

20 . 000 

-2.29 

5.11 

-.0450 

-.23 

2.68 

5.22 

3.35 

1.04 

322. 

21 .COO 

-2.76 

6.25 

-.0027 

.03 

2.47 

6.04 

4.03 

.96 

309. 

22.000 

-3.52 

7.28 

-.0759 

-.03 

2.21 

7.16 

4.35 

.98 

306. 

23.000 

-4.24 

8.66 

-.1078 

-.  18 

2.22 

6.64 

4.81 

.63 

300. 

24. COO 

-4.98 

9.7C 

-.1236 

-.18 

2.29 

9.88 

5.13 

.28 

303. 

25.000 

-5.  13 

10.71 

-.0616 

-.08 

2  32 

10.74 

5.55 

.31 

297. 

26.000 

-5.22 

12,04 

.024  3 

-.13 

2.06 

1 1  .69 

6. GO 

.29 

203. 

27.000 

-4.93 

13.49 

.0509 

.01 

3.43 

12.99 

6.98 

.16 

257. 

2d .  r.ao 

-4.25 

14.79 

.1183 

.  15 

3.95 

14.42 

6.60 

.03 

242. 

29.020 

-3.02 

15.60 

.  1  334 

.25 

4 . 24 

14.95 

6.80 

.  17 

218. 

30 . 000 

-2.64 

16.37 

.0781 

.35 

4.3* 

15.23 

7.74 

•  3-< 

172. 
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KAPCH 


note 

STATION  • 

2 

KM 
.til 
1 .000 
3.000 
3.000 
4.000 
5.000 
6.000 
7.000 
8.000 
9.000 
10.000 
1 1 .000 
13.000 
13. COO 
1<*.000 
IS. COO 
Ifi.COO 
17.000 
10.000 
19.000 
30 . 000 
31.000 

3; .  ooo 
33.000 
3*t .  000 
30.000 
36.000 
37 . 000 
30 . 000 
39.000 
30 . 000 


TA91X 

STATION 

2 

KM 

.lit 

1.000 

3.000 

3.000 

<4.000 

5.000 

6.000 

7.000 

8.000 

9.000 

10.000 

11.000 

13.000 

13.000 

1*4.000 

15.000 

16.000 

17.000 

18.000 

19.000 

30.  COO 

31.  GOO 
33.000 
33.000 
3H.000 
35.000 
36.000 
37.000 
38.000 
39.000 
30.000 
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TAGUAC 

MEAN  U 

S.O.  U 

RIU.V) 

TCAN  V 

M/S 

M/S 

M/S 

-*4.78 

1  .89 

.0367 

-1.33 

-9.71 

3.  *45 

.0958 

-1.59 

-7.13 

3.83 

.1793 

-1.09 

-5.55 

3.96 

.3169 

-.75 

-*4.35 

3.88 

.1331 

-.71 

-3.7*4 

*4.33 

.0937 

-.71 

-3.31 

*4.76 

.0870 

-.59 

-3.  *48 

5.15 

-.0609 

-1 .05 

-1.36 

5.98 

- .3039 

-1 .07 

.0*4 

5.75 

-.1913 

-.65 

1.13 

5.99 

-.1731 

-.07 

1.97 

6.93 

-.1595 

1.09 

3.60 

6.58 

-.1197 

3.79 

3. >41 

6.58 

-.1531 

3.99 

1.71 

6.73 

-.1993 

3.96 

-.06 

6.93 

-.3993 

3.1! 

-3.99 

6.93 

-.9989 

3.76 

-3.95 

6.36 

-.3616 

1.76 

-3. 18 

9.59 

-.1396 

.99 

-1.87 

9.33 

-.oeeo 

-.30 

-1.17 

9.79 

-.0930 

-.31 

-1.65 

5.79 

-.0651 

.07 

-3.90 

6.96 

-.0893 

.  15 

-3.91 

8.33 

-.0691 

-.  to 

-*4.57 

9.33 

-.1333 

-.33 

-5.38 

10.03 

-.1035 

-.31 

-5.73 

11.07 

-.0563 

-.93 

-5.95 

13.11 

-.0937 

-.93 

-5.6*4 

13.79 

-.0313 

-.93 

6.03 

13.09 

.  0006 

-.35 

-*j.7S 

13.93 

-.01  lb 

•  GO 

1.  9 

Ml NO  STATISTICAL  PARAMETERS. 
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TAGUAC 

TCAN  U 

S.O.  U 

RIU.V) 

TCAN  V 

M/S 

M/S 

M/S 

-5.98 

1.89 

.0933 

-.96 

-10.63 

3.95 

.0509 

-.51 

-8.61 

3.36 

.1333 

-.63 

-6.90 

3.80 

.3016 

-.70 

-6.13 

370 

.3195 

-.90 

-5.98 

3.88 

.1899 

-1.35 

-5.38 

9.87 

.3399 

-1.33 

-9. 13 

5.61 

.1099 

-1.99 

-1.88 

6.95 

.0670 

-1 .38 

1.93 

7.05 

.1573 

-1.33 

9.07 

7.99 

.1778 

-.93 

6.09 

8.03 

.1376 

.  10 

7.18 

8.53 

.1167 

t.39 

7.36 

8.73 

.06-49 

1.39 

6.99 

8.39 

.0096 

.37 

9.03 

7.59 

-.0659 

-1.03 

.03 

<-..31 

-.1108 

-1.36 

-3.71 

9.93 

-.1377 

-1 .99 

-9.80 

9.39 

-.0591 

-1.33 

-9.93 

9.11 

-.0397 

-.93 

-3.79 

9.99 

.0179 

-.53 

-9.06 

6.35 

.0708 

-.35 

-5.91 

7.99 

-.0703 

-.13 

-6.79 

8.53 

-.0915 

-.08 

-7.89 

9.18 

.0335 

-.06 

-9.16 

9.65 

.1038 

-.07 

-10.19 

10.37 

.0807 

.36 

-10.95 

10.80 

.0699 

.39 

-11.53 

10.61 

.0910 

.36 

-11.88 

10.56 

.0999 

.59 

-13.79 

11.08 

.1973 

.99 

S.O.  V 

MEAN  MS 

S.O.  MS 

SKEM  MS 

M/S 

M/S 

M/S 

1.69 

5.33 

1.78 

.31 

3.  18 

10.93 

3.30 

.37 

3. 89 

7.98 

3.30 

.10 

3.36 

6.90 

3. 17 

.93 

3.67 

6.17 

3.13 

.61 

9.08 

6.35 

3.  16 

.58 

9.68 

6.66 

3.38 

.67 

9.78 

6.80 

3.31 

.58 

9.66 

6.63 

3.37 

.57 

9.91 

6.30 

3.63 

.71 

9.90 

6.93 

3.81 

.71 

9.83 

7.10 

9.33 

.73 

5.  36 

8.09 

9.55 

.83 

5.63 

8.53 

9.85 

1.05 

5.79 

8.95 

5.03 

1.11 

5.68 

7.97 

5.  13 

.98 

5.93 

8.30 

5.00 

1.16 

9.37 

7.70 

9.39 

1.03 

3.81 

5.99 

3.93 

.56 

3.99 

9.63 

3.93 

.50 

3.97 

9.95 

3.99 

.93 

3.57 

5.53 

3.51 

.96 

3.30 

6.51 

9.39 

.85 

3.18 

7.97 

5.08 

.80 

3.96 

9.16 

5.30 

.55 

3.90 

10.05 

5.83 

.96 

3.76 

10.86 

6.70 

.91 

3.35 

11.69 

7.97 

.91 

3.73 

IP.  33 

7.57 

.96 

9 

13.00 

7.91 

.93 

9.96 

13.11 

6.93 

.57 

APRIL 


S.O.  V 

MEAN  MS 

S.O.  MS 

SKEM  MS 

M/S 

M/S 

M/S 

1.33 

5.76 

1.73 

.05 

3.36 

10.93 

3.70 

-.56 

3.91 

9.05 

3.10 

-.03 

3.79 

7.73 

3.31 

.02 

3.83 

7.16 

3.96 

.36 

3.88 

7.07 

3.27 

.99 

3.98 

7.39 

3.56 

.39 

3.93 

7.35 

3.67 

.69 

9.31 

7.33 

3.66 

.65 

9.79 

7.57 

9.30 

.67 

5.09 

8.69 

9.96 

.63 

6.17 

10.35 

5.63 

.55 

6.91 

11.67 

6.12 

.50 

6.53 

11.65 

6.12 

.59 

5.60 

10.98 

6.79 

.63 

9.06 

8.68 

9.79 

.93 

3.95 

6.96 

3.77 

.78 

3.51 

5.97 

3.23 

.81 

3.  13 

6.31 

3.56 

.89 

3.96 

5.83 

3.05 

.61 

3.33 

5.97 

2.96 

.97 

3.  16 

6.71 

3.89 

.67 

3.08 

8.39 

9.66 

.63 

3.35 

9.60 

5.62 

.90 

3.37 

10.52 

6.36 

.31 

3.97 

11.33 

7.91 

.35 

3.70 

13.31 

8.  15 

.36 

3.19 

13.17 

8.53 

.31 

3.39 

13.90 

7.99 

.39 

3.55 

19.95 

7.52 

.36 

3.83 

15.58 

7.62 

.31 

MAY 


TABLE  t.  5  W1N0  STATISTICAL  PARA1CTERS. 


STATION 

•  912170 

TAOUAC 

z 

MEAN  U 

S.D.  U 

R1U.VI 

1C  AN  V 

S.D.  V 

MEAN  WS 

S.D.  WS 

SKEW  WS 

NOBS 

KM 

M/S 

M/S 

M/S 

M/S 

M/S 

M/S 

•  111 

-4.11 

2.44 

.3371 

.00 

1 .87 

4.76 

1.91 

.20 

368. 

1.000 

-7.08 

4.63 

.2992 

.43 

2.95 

8.18 

3.66 

.42 

415. 

a.  ooo 

-5.36 

4.  12 

,2162 

.46 

2.72 

6.54 

3.26 

.89 

416. 

3.000 

-4.49 

4.31 

.1133 

.24 

2.99 

6.16 

3. 12 

.89 

414. 

4.000 

-3.83 

3.66 

.1016 

-.23 

2.94 

5.40 

3.02 

1.07 

411  . 

5.000 

-3.36 

4.07 

.  1662 

-.49 

2.90 

5.22 

3.04 

1 . 17 

411. 

6.000 

-3.05 

4.11 

.1410 

-.71 

3.24 

5.30 

3.01 

1 .02 

408. 

7.CC0 

-2.11 

4.69 

-.0760 

-.98 

3.41 

5.44 

3.07 

.66 

403. 

8.000 

-.65 

5.66 

- .0050 

-1.62 

3.48 

5.89 

3.36 

.66 

4C6. 

9.  COO 

1.49 

6.87 

.  ICO* 

-1.89 

4.29 

7.31 

4.21 

.77 

406. 

10.000 

3.53 

8.60 

.1956 

-1  .58 

4  .  'h 

9.16 

5.24 

.79 

405. 

11 .000 

5.34 

9.93 

.2359 

-1.07 

6.33 

11.37 

6.34 

.68 

403. 

12. COO 

6.14 

11.39 

.3496 

-.74 

7.13 

13.05 

6.95 

.54 

401  . 

13.000 

6. 15 

12.04 

.  3472 

-1.33 

7.07 

13.56 

7.08 

.53 

399. 

14.000 

5.17 

11.50 

.3208 

-2.08 

6.41 

12.60 

6.75 

.66 

399. 

15.000 

2.37 

9.81 

.2610 

-2.91 

5  21 

10.35 

5.49 

.59 

39G. 

16.000 

-1.48 

7.54 

.2384 

-3.50 

3.82 

8.33 

4.05 

.48 

388. 

17.000 

-4.19 

6.20 

.0770 

-3.21 

3.08 

7.90 

3.66 

.37 

346. 

18.000 

-6.05 

4.48 

-.0568 

-2.21 

2.56 

7.57 

3.30 

.15 

348. 

19.000 

-7.01 

3.82 

-.0477 

-1.25 

2.36 

7.74 

3.30 

.  18 

347. 

20.000 

-8. 14 

4.29 

-.0143 

-.34 

2.06 

8.56 

3.95 

.40 

343. 

21.000 

-9.36 

5.26 

.0209 

-.21 

2.06 

9.78 

4.88 

.28 

335. 

22.000 

-10.31 

5.69 

-.054! 

-.11 

2,22 

11.25 

5.47 

.11 

333. 

23.000 

-12.41 

6.43 

-.0060 

-.09 

2.40 

12.76 

6.21 

.09 

325 . 

24.000 

-13.67 

7.35 

-.0046 

-.14 

2.55 

14.08 

7.02 

.08 

321. 

20.000 

-15.44 

7.69 

-.0919 

-.08 

2.31 

15.69 

7.53 

.23 

316. 

26.000 

-16.93 

7.91 

-.1056 

-.37 

2.59 

17.19 

7.77 

.13 

307. 

27.000 

-10.31 

7.91 

-.0623 

-.22 

2.74 

18.58 

7.76 

.07 

209. 

20 . 000 

-18.87 

8.26 

-.0665 

.00 

3.02 

9.19 

8.07 

.06 

265. 

29.000 

-18.77 

8.45 

-.0806 

.19 

3.43 

!  9 . 1 9 

8.20 

.22 

237. 

30 . 000 

-10.74 

B.’iS 

-.0475 

.  10 

3.21 

19.05 

8.32 

.31 

!  79. 

TABLE 

1.  6 

WIND  ST AT! 

IST1CAL  PARAfCTERS. 

JCfC 

STATION 

•  912170 

TAGUAC 

Z 

MEAN  U 

S.O.  U 

R1U.V1 

fCAN  V 

S.D.  V 

MEAN  WS 

S.D.  WS 

SKEW  WS 

Noes 

KM 

M/S 

M/S 

M/S 

M/S 

M/S 

M/S 

.111 

-4.02 

1.89 

.1266 

.23 

1.43 

4.37 

1.64 

.11 

378. 

1.000 

-7.57 

3.49 

.0713 

1.18 

2.59 

8.34 

2.83 

-.03 

411. 

2.000 

-6.33 

3.06 

.1299 

.65 

2.55 

7.08 

2.49 

.32 

413. 

3.000 

-6.40 

3. 12 

.1272 

.69 

2.73 

7.16 

2.70 

.40 

413. 

4.000 

-5.86 

3.08 

.0899 

.  15 

2.88 

6.64 

2.84 

.24 

413. 

5.U00 

-5.53 

3.28 

.1151 

-.11 

2.78 

6.41 

2.82 

.22 

413. 

6.000 

-5.20 

3.58 

.1133 

-.26 

3.20 

6.42 

2.99 

.22 

412. 

7 . 000 

-4.36 

3.93 

.1697 

-.22 

3.51 

6.22 

2.85 

.27 

41  1  . 

8.000 

-3.39 

4.80 

.1619 

-.27 

4.20 

6.46 

3.23 

.50 

407. 

9. COO 

-2.07 

6.00 

.2556 

-.50 

5.23 

7.20 

3.99 

.62 

406. 

10.000 

-.79 

6.94 

.2842 

-.38 

6.29 

8.23 

4.53 

.72 

406. 

11.000 

.16 

8.23 

.2923 

-.23 

7.29 

9.68 

5. 19 

.52 

4C5. 

12.000 

1.17 

9.67 

.2914 

-.75 

8.37 

11.38 

5.97 

.46 

403. 

13.000 

1.57 

10.70 

•  2B37 

-1.36 

9.  19 

12.67 

6.52 

.46 

403. 

14.000 

1.43 

10.92 

.2871 

-2.01 

9.38 

12.79 

7.01 

.60 

403. 

15.000 

-.30 

10.30 

.3357 

-2.68 

7.24 

11.13 

6.46 

.87 

3S3. 

16.000 

-4.35 

7.66 

.2882 

-2.59 

4.68 

8.74 

5.49 

1.23 

377. 

17.000 

-7.20 

5.41 

.2298 

-2.29 

3.  10 

8.60 

4.69 

.89 

357. 

18.000 

-9.83 

3.68 

.1281 

-1.63 

2.24 

10.25 

3.56 

.25 

358. 

19.000 

-11.43 

3.73 

-.0262 

-1.22 

1  .93 

11.68 

3.66 

.38 

356 . 

20 . 000 

-12.91 

4.02 

.1113 

-.56 

2.05 

13.10 

3.98 

.41 

350. 

21.000 

-14.53 

4.89 

.0181 

-.20 

2.  18 

14.71 

4.84 

.26 

343. 

22.000 

-16.22 

5.72 

-.0575 

-.08 

2.07 

16.37 

5.67 

.27 

350. 

23.000 

-17.66 

6.45 

-.0626 

-.17 

2.35 

17.83 

6.40 

.19 

324. 

24.000 

-19. 12 

6.77 

-.1109 

-.24 

2.57 

19.31 

6.72 

.07 

318. 

25.000 

-20.48 

7.11 

-.0694 

-.02 

2.26 

20.62 

7.05 

.00 

323. 

26.000 

-21.74 

7.38 

-.0861 

-.08 

2.43 

21  .90 

7.30 

-.10 

308. 

27.000 

-23.09 

7.34 

-.1027 

.03 

2.59 

23.26 

7.26 

-.27 

287. 

28.000 

-24. 10 

7.67 

-.0168 

.04 

2.97 

24.30 

7.60 

-.16 

264. 

29.000 

-24.90 

7.70 

-.0294 

-.06 

3.03 

25.10 

7.67 

-.07 

232. 

30. COO 

-25.53 

7.63 

-  .0839 

-.06 

3.45 

25.79 

7.52 

-.03 

178. 

59 
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TABLi 

I.  7 

MIND  STAT 1 

:ST!CAL  PARAfCTE'TS, 

Jl 

ILY 

STATION 

•  912170 

TAOUAC 

Z 

MEAN  U 

S.O.  U 

R(U.V) 

tCAN  V 

s.o.  v 

MEAN  MS 

S.O.  MS 

SKEM  MS 

nobs 

KM 

M/S 

M/S 

M/S 

M/S 

M/S 

M/S 

.111 

-2.32 

2.32 

.0909 

.57 

l  .74 

3.38 

1.63 

.54 

409. 

1.000 

-4.44 

5.04 

.0968 

1.61 

3.42 

6.76 

3.69 

.80 

423. 

2.000 

-4.02 

4.54 

.2503 

I. 37 

1.20 

6. 19 

3.23 

.82 

425. 

3.000 

-4.43 

4.53 

.2333 

1.29 

2.35 

6.53 

3.25 

.48 

425. 

4.000 

-4.22 

4.64 

.3016 

1 .34 

3.60 

6.58 

3.27 

.61 

427. 

5.000 

-4.  16 

4.53 

.3254 

1 .  1 1 

3.92 

6.64 

3.30 

.58 

427. 

G.000 

-4.00 

4.51 

.2925 

.78 

3.88 

6.36 

3.43 

.74 

42*. 

7.000 

-3.81 

4.22 

.2885 

.21 

3.71 

5.96 

3.26 

.67 

421  . 

8.000 

-3.54 

4.25 

.2302 

-.25 

3.67 

5.93 

3.23 

.59 

421  . 

9.0C0 

-3.34 

4.86 

.1717 

-.61 

4.23 

6.39 

3.50 

.83 

421  . 

10.000 

-3.17 

5.44 

.1650 

-1.05 

4.87 

6.96 

3.99 

1.17 

419. 

11.000 

-3.22 

6.53 

.  16*4 

-1.66 

5.68 

8.11 

4.77 

1 .21 

415. 

12. GOO 

-3.30 

7.90 

.1787 

-2.95 

6.79 

9.88 

5.60 

.91 

415. 

13.000 

-3.83 

8  84 

.1646 

-4.29 

7.84 

1  1  .44 

6.43 

.99 

409. 

14.000 

-4.71 

9.42 

.1341 

-5.23 

8.35 

12.60 

7.01 

.74 

405. 

15.000 

-6.70 

8.69 

.2663 

-4.99 

7.  14 

12.21 

6.85 

.65 

396. 

16.000 

-9.27 

6.42 

.2897 

-3.74 

4.42 

1 1  .40 

5.53 

.41 

378. 

17.000 

-1 1 .59 

4.25 

.2537 

-2.60 

2.67 

12.20 

4.20 

-.16 

358. 

18.000 

-14.05 

3. 15 

.0734 

-1.43 

2.26 

14.31 

3.13 

.26 

359. 

19.000 

-15.71 

3.42 

.0337 

-1 .01 

2.20 

15.90 

3.41 

.34 

353. 

20 . 000 

-16.92 

3.88 

-.0077 

-.39 

2.11 

17.06 

3.84 

.  16 

346. 

21.000 

-18. 15 

4.82 

.0456 

-.24 

2.45 

18.33 

4.76 

.26 

336. 

22.000 

-19.48 

5.78 

-.0058 

-.07 

2.22 

19.63 

5.68 

.22 

313. 

23.000 

-20. 8G 

6.65 

-.0268 

-.02 

2.41 

21.03 

6.55 

.18 

303. 

24 . 000 

-22.15 

7.  14 

-.0839 

-.23 

2.57 

22.32 

7.08 

.09 

295. 

25.000 

-23.99 

7.14 

-.0712 

-.09 

2.  19 

24.09 

7. 12 

-.13 

298. 

26.000 

-25.31 

7.  17 

-.0499 

.01 

2.38 

25.43 

7.  14 

-.30 

288 

27.000 

-26.72 

7.00 

-.0883 

.02 

2.77 

26.87 

6.94 

-.15 

266. 

28.003 

-27.96 

7.32 

-.0878 

.16 

3.14 

28.15 

7.27 

-.  12 

246. 

29.000 

-29.38 

7.79 

-.0372 

.17 

2.83 

29.54 

7.73 

-.  14 

2)0. 

30.000 

-30.18 

7.28 

-.0657 

.41 

3.08 

30.35 

7.23 

-.21 

145. 

TABLE  1 .  8 

STATION  -  912170 
Z  KEAN  U 


MIND  STATISTICAL  PARAMETERS, 
TAGUAC 


KM 
.III 
1. 000 
2.000 
3.000 
4.000 
5.000 
6.000 
7.000 
8.000 
9.000 
10.000 
11.000 
12.000 
13.000 
14.000 
15.000 
16.000 
17.000 
18.000 
19.000 
20.000 
21 .000 
22.000 
23.000 
.000 
.000 
.000 
.000 
.000 
.000 
30.0J0 


M/S 
-I. HO 
-2.08 
-2.18 
-2.85 
-3.06 
-3.03 
-3.01 
-2.87 
-2.66 
-2.59 
-2.64 
-2.66 
-2.87 
-3.02 
-3.72 
-5.90 
-8.94 
-11.50 
-14.29 
-15.80 
-16.99 
-18.71 
-19.98 
-21.09 
-22.21 
-23.78 
-25. 14 
-26.70 
-27.72 
-28.59 
-30.29 


s.o.  u 

M/S 

2.81 

6.71 
6.38 
5.93 
5.36 

5.43 

5.29 
5.20 
5.33 
5.66 
6.25 

7.30 
8.17 
9.05 

9.84 

9.77 
6.98 

4.77 

3.82 
3.80 
4.41 
5.09 
5.93 
6.48 

6.72 
6.69 

6.84 

6.86 

7.43 
7.65 

7.82 


RCU.V1 

.2512 

.3023 

.3228 

.2657 

.2616 

.2671 

.1971 

.1933 

.2211 

.2351 

.2686 

.3018 

.2994 

.2964 

.3106 

.2902 

.2070 

.2633 

.0601 

.0760 

.0150 

-.0544 

-.0426 

-.0232 

-.1371 

-.0738 

-.1852 

-.1502 

.0552 

.0074 

.0042 


MEAN  V 
M/S 
.80 
2.13 
1.82 
1.77 
1.75 
1.58 
1.18 
.72 
.39 
-.19 
-.98 
-1.79 
-3.28 
-4.90 
-6.  13 
-5.93 
-3.96 
-2.58 
-1 .27 
-.76 
-.53 
-.23 
■  00 
.07 
.09 
-.06 
-.11 
.08 
.27 
.37 
.  16 


AUOUST 


S.O.  V  MEAN  US  S.O.  US 
M/S  M/S  M/S 

1.98  3.35  1.78 


4.41  7 

4.29  6 

4.28  6 

4.54  6 

4.80  6 

4.67  6 

4.13  6 

4.24  6 

4.59  6 

4.93  7 

5.94  8 

7.31  10 

8.80  12 

9.56  13 

8.24  13 

4.91  II 

2.87  12 

2.22  14 

2.25  15 

2.28  17 

2  22  18 

2.30  20 

2.47  21 

2.63  22 

2.36  23 

2.57  25 

2.89  26 

3.01  27 

3.12  28 

3.31  30 


.23  4.57 
■8B  4.45 
.96  4.04 
.84  3.85 
.88  4.09 
.66  3.99 
.30  3.61 
.37  3.59 
.81  3.66 
.53  3.82 
.73  4.75 
.23  5.86 
.04  6.87 
.61  7.36 
.49  7.14 
.41  6.17 
.15  4.71 
.53  3.79 
.97  3.79 
.15  4.39 
.85  5.06 
.13  5.88 
.27  6.39 
.39  6.63 
.91  6.64 
.29  6.79 
.87  *.79 
.91  7.32 
.79  7.52 
>9  7.76 


SKEW  MS 

.95 
1.65 
1.91 
1.34 
1.54 
1.64 
1 .26 
.93 
.99 
.72 
.50 
.82 
.89 
.88 
.74 
.59 
.46 
.41 
.22 
.31 
.50 
.27 
.32 
.  1 1 
-.03 
-.03 
-.05 
-.16 
-.49 
-.60 
-.22 


NOSS 


398. 

422. 

425. 

424. 

425. 
425. 
425. 

422. 

423. 
422. 
421  . 
418. 
417. 
417. 
415. 
409. 
392. 
360. 
359. 
352. 
345. 
334. 
318. 
309. 
304. 
2G9. 
277. 
254. 
229. 
196. 
140. 
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table 

l .  9 

MIND  j.aY' 

1ST1CAL  PARAMETERS 

station 

•  91317C 

TAOUAC 

2 

ME  A,*1  ij 

s.o.  u 

RfU.V) 

MEAN  V 

KM 

M/S 

M/S 

M/S 

•  Ill 

-1.47 

3.03 

.1470 

.30 

1 .000 

-3.79 

6.35 

.1788 

1.04 

3.000 

-3.03 

5.77 

.3384 

.90 

3.000 

-::.78 

5.57 

.31  17 

.09 

■4.000 

-3.95 

5.68 

.2381 

.89 

5.000 

-3.91 

5.70 

.3985 

.00 

6.000 

-3  96 

5.54 

.  3450 

.63 

7.000 

-3.57 

5. 15 

.3392 

.37 

a  coo 

-3.35 

5.31 

.2391 

.  12 

9.000 

-2. 7*4 

5.37 

.1563 

-.16 

10.000 

-3.33 

5.63 

.1557 

-.53 

1! .000 

-3.00 

6.45 

.  1  167 

-.77 

13. COO 

-1.93 

7.40 

.0248 

-1.33 

1 f . 000 

-1  .81 

8749 

-.0249 

-1.86 

1  1  . ooo 

-3.61 

9.08 

.0170 

-3.43 

15. COO 

-5.09 

8.68 

.1344 

-3.74 

16. COO 

-8.3-4 

6.51 

.1489 

-2.36 

17.1:00 

-10. >43 

•4 . 29 

.0698 

-1  .46 

19.  ino 

-13.08 

3.37 

.0413 

-.69 

19.000 

-13.32 

3.72 

.0677 

-.47 

30 . C CO 

-  1 H . 44 

4.42 

.0368 

-.30 

=•  1  . 1  C’O 

-15.98 

5  14 

-.0023 

-.03 

r  ' .  0 0 

-17.68 

5.07 

-.0737 

.33 

33.000 

-19.03 

6.  18 

-.0761 

.38 

34.000 

-30.33 

6.47 

-.0903 

.  13 

35 . 000 

-21.96 

7.03 

.0344 

.02 

■TO.  COO 

-33.91 

7.43 

.03/9 

.00 

27.000 

-39 . 56 

8.06 

.0333 

-  .08 

36  000 

-25.88 

8.31 

.0733 

.  31 

39  000 

-25.57 

8.44 

.1139 

.78 

30.000 

-27.93 

b  e» 

.  0 1  34 

.46 

TA3LE  t.  10  MIND  STATISTICAL  PARAfCTERS. 
STATION  -  913170  TAOUAC 


S.O.  V 

MEAN  MS 

S.O.  MS 

SKEW  MS 

NOBS 

M/S 

M/S 

M/S 

1.87 

3.36 

1.78 

.59 

389. 

3.98 

6.90 

4.01 

1.31 

407. 

3.64 

6.41 

3.91 

1.71 

400. 

3.60 

6.79 

3.68 

1 .03 

408. 

3  89 

7.03 

3.79 

.95 

409. 

3.93 

6.96 

3.91 

1 .34 

409. 

4.03 

6.91 

3.90 

1.05 

408. 

4.00 

6.50 

3.61 

1.15 

408. 

4.10 

6.46 

3.66 

1 . 10 

408. 

4.83 

6.63 

3.81 

1 . 14 

406. 

5.30 

7.04 

3.84 

.98 

406. 

6.13 

8.00 

4.41 

.83 

405. 

7.13 

9.33 

4.88 

.83 

404. 

7.89 

10. 60 

5.34 

.65 

393. 

8.  13 

11.37 

5.83 

.71 

397. 

6.93 

10.69 

6.49 

1.35 

390. 

4.  15 

9.98 

5.88 

1.06 

364. 

3.70 

10.90 

4.18 

.50 

335. 

3.36 

12  30 

3.34 

.13 

334. 

3.06 

13.50 

3.69 

.35 

339. 

3.06 

14.60 

4.38 

.  18 

324. 

2  23 

16. 14 

5.13 

.07 

316. 

2.13 

17.83 

5.83 

.13 

309. 

3.31 

19.  19 

6.13 

.07 

303. 

3.57 

30.43 

6.38 

.11 

397. 

3.19 

22.08 

6.99 

.  10 

391. 

3.64 

33.10 

7.33 

-.33 

375. 

3.83 

34.75 

7.96 

-.40 

357. 

3.83 

26.08 

8.06 

-.60 

338. 

3.9 

36.78 

8.31 

-.60 

316. 

3.65 

28.31 

8  71 

-.55 

145. 

OCTOBER 


Z 

MEAN  U 

S.D.  U 

R(U.V) 

MEAN  V 

KM 

M/S 

M/S 

M/S 

.III 

-3.94 

3.78 

.1680 

.33 

1.000 

-6.00 

6.34 

.0698 

1 .36 

3.000 

-5.48 

5.68 

.0633 

1.03 

3.000 

-5.89 

5.71 

.0737 

1.08 

4.000 

-5.89 

5.79 

.0655 

.76 

5.000 

-5.91 

5.93 

.0570 

.55 

6.000 

-5.78 

5.90 

.0797 

.48 

7.000 

-5.05 

5.91 

.1575 

.31 

8.000 

-4. 19 

6.03 

.1740 

.  13 

9.000 

-3.88 

6.39 

.1731 

.33 

10.000 

-1.76 

6.89 

.1685 

.53 

11.000 

-.87 

8.C6 

.1906 

.40 

12.000 

.04 

8.83 

.1530 

.09 

13.000 

.03 

9.41 

.0987 

-.61 

14.000 

-1.13 

9.  16 

.0963 

-1.56 

15.000 

-4.61 

8.49 

.1118 

-3.68 

IS. 000 

-8.36 

7.07 

-.0103 

-3.56 

17.000 

-9.89 

5.31 

.0313 

-3.17 

18.000 

-10.35 

4.36 

•  CC06 

-.98 

19.000 

-10.38 

4.55 

-.1178 

-.44 

30.000 

-11.61 

5.  14 

-.0969 

-.13 

21 .000 

-13.08 

6.09 

-.0947 

.15 

23.000 

-14.48 

6.63 

.0159 

.33 

33.000 

-15.63 

7.35 

-.1373 

.01 

34.000 

-16.69 

8.11 

-.1095 

-.09 

35.000 

-17.01 

8.77 

-.0855 

-.03 

36.000 

-18.63 

9.40 

-.1153 

.13 

37.000 

-19.36 

9.55 

-.3064 

.37 

38.000 

-18.87 

9.65 

- .3500 

.71 

39.000 

-17.77 

9.73 

- .  1 3S8 

.76 

30 . COO 

-17.39 

10.55 

-.0751 

.43 

S.D.  V 

MEAN  MS 

M/S 

M/S 

1.98 

4.05 

4.34 

8.73 

3.85 

7.89 

4.  19 

8.35 

4.37 

8.34 

4.58 

8.40 

4.56 

8.37 

4.67 

7.89 

4.64 

7.56 

5.34 

7.68 

6.11 

8.34 

7.38 

9.54 

8.11 

10.53 

8.44 

11.05 

8.33 

10.71 

6.77 

10. 17 

5.30 

10.66 

3.64 

tO.  83 

3.44 

10.74 

3.36 

10.73 

3.  15 

11.86 

3.03 

13.30 

3.39 

14.76 

3.07 

15.89 

3.51 

17.05 

3.  14 

18.09 

3.68 

19.00 

3.03 

19.67 

3.30 

19.41 

3.77 

18.51 

3.60 

19.21 

D.  MS 

SKEW  MS 

M/S 

1.98 

3.00 

4.39 

.83 

3.99 

1.07 

4.31 

.95 

4.47 

.88 

4.55 

.89 

4.38 

.78 

4.46 

.95 

4.36 

.83 

4.30 

.97 

4.30 

.81 

5.38 

.79 

5.73 

.90 

6.15 

.99 

6.33 

1.03 

6.53 

1.16 

6.33 

.91 

5. 17 

.66 

4.15 

.30 

4.34 

.  13 

5.03 

.30 

5.95 

.  18 

6.40 

.07 

6.96 

.03 

7.73 

.16 

8.44 

.  14 

9.01 

-.01 

9. 19 

-.30 

9. 17 

-.08 

9.07 

.08 

9.76 

.33 

NOBS 

380. 

>417. 

418. 

ma. 

419. 

ma. 

ma. 

•♦17. 

ms. 

m3. 

mo. 

mo. 

*409. 

>ioe. 

>406. 

■400. 

369. 

330. 

318. 

310. 

303. 

393. 

300. 

379. 

377. 

366. 

35m 

336. 

313. 

174. 

134. 


TABLE  !.  13  WIND  STATISTICAL  PARAMETERS,  AWUAL 

STATION  •  912170  TAGUAC 


z 

1C  AN  U 

S.D.  U 

R(U.V) 

1C  AN  V 

S.D.  V 

MEAN  US 

S.D.  US 

SKEW  US 

NOBS 

KM 

M/S 

M/S 

M/S 

M/S 

M/S 

M/S 

•  111 

-3.66 

2.65 

.2808 

-.30 

1 .92 

4.51 

1.97 

.34 

4514. 

1.000 

-7.34 

5.70 

.2779 

.  17 

3.91 

9.21 

4. 10 

.41 

4960. 

2.000 

-6. 1 1 

5.23 

.2526 

.27 

3.61 

7.91 

3.9U 

.81 

4976. 

3.000 

-5.67 

5.09 

.1967 

.38 

3.72 

7  .‘'9 

3.80 

.75 

4973. 

4.000 

-5.27 

4.96 

.1678 

.35 

3.86 

7.27 

3.82 

.85 

4974. 

5.000 

-5.16 

5.10 

.1228 

.23 

4.09 

7.33 

3.96 

.90 

4975. 

6.000 

-4.99 

5.34 

.0719 

.20 

4.37 

7.45 

4.13 

.84 

4961. 

7.000 

-4.36 

5.43 

.0238 

-.12 

4.42 

7.22 

3.99 

.87 

4944  . 

B.000 

-3.42 

5.77 

-.0193 

-.39 

4.51 

7.07 

3.94 

.89 

4916. 

9.000 

-2.12 

6.37 

.0169 

-.55 

4.82 

7.19 

4.12 

.92 

4907. 

10.000 

-.99 

7.06 

.0842 

-.47 

5.18 

7.68 

4.36 

.91 

4887. 

11.000 

-.09 

8.01 

.1518 

-.16 

6.03 

8.70 

4.99 

.91 

4870. 

12.000 

.44 

8.76 

.1862 

.04 

7.01 

9.76 

5.55 

.92 

4849. 

13.000 

.39 

9.28 

.1670 

-.08 

7.78 

10.54 

5.98 

.94 

4819. 

14.000 

-.37 

9.37 

.1414 

-.41 

8.13 

10.74 

6.24 

.93 

4801  . 

15.000 

-2.61 

8.91 

.1035 

-.76 

7.39 

10.21 

6. 10 

.95 

4740. 

16.000 

-5.82 

7.44 

-.0123 

-.58 

5.77 

9.53 

5.65 

.90 

4569. 

17.000 

-7.72 

6.17 

-.0336 

-.48 

4.39 

9.62 

4.98 

.58 

4186. 

18.000 

-8.69 

5.61 

.0439 

-.47 

S.9* 

9.64 

4.80 

.35 

4193. 

19.000 

-8.65 

6.23 

.0431 

-.48 

2.39 

9.56 

5.31 

.42 

4135. 

20.000 

-0.79 

7.30 

.0169 

-.24 

2.26 

10.04 

5.92 

.49 

4077. 

21.000 

-9.68 

8.43 

.0146 

-.04 

2.27 

11.16 

6.73 

.43 

3940. 

22.000 

-10.84 

9.14 

-.0316 

.07 

2.18 

12.37 

7.27 

.43 

3857. 

23.000 

-11.89 

9.95 

-.0465 

.01 

2.31 

13.65 

7.72 

.40 

3770. 

24.000 

-12.78 

10.64 

-.0413 

-.05 

2.50 

14.73 

8.13 

.35 

3746. 

25.000 

-13.71 

11.48 

-.0375 

-.04 

2.35 

15.85 

8.62 

.30 

3676. 

26.0C0 

-14.26 

12.36 

-.0239 

-.02 

2.70 

16.73 

9. 13 

.20 

3539. 

27.000 

-14.66 

13.49 

-.0178 

.04 

3.01 

17.71 

9.61 

.13 

3269. 

23.000 

-14.48 

14.57 

.0002 

.17 

3.36 

18.37 

9.80 

.19 

3028. 

23.000 

-14.10 

15.48 

.0232 

.29 

3.68 

18.80 

9.90 

.28 

2604. 

37.000 

-13.56 

16.60 

.0497 

.32 

3.80 

19.36 

10.  10 

.32 

1099. 

TtCWICOYNAMC  STATISTICAL  PARAJCTERS 
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TABLE  IV. t  HYDROSTATIC  MODEL  ATMOSPHERE.  TABLE  IV.  2  HYDROSTATIC  MODEL  ATMOSPHERE. 

STATION  •  912170  TAOUAC  STATION  -  912170  TAGUAC 


Z 

GCO.  HT. 

P 

0 

TV 

Z 

GEO.  HT. 

P 

D 

TV 

KM 

KM 

MS 

G/M3 

DEO  K 

KM 

KM 

M8 

G/M3 

OEO  K 

.000 

.000 

1012.5000 

1169.0000 

301.73 

.000 

.000 

1013.1000 

1 170.0000 

301 .65 

.111 

.111 

999.9000 

1 157.0000 

301 .00 

.111 

.1)1 

1000.5000 

1158.0000 

300.9*7 

1  .000 

.997 

903.2100 

1068.0000 

294.64 

1.000 

.997 

903.6700 

1070.0000 

294.26 

2.000 

1 .995 

803.7700 

967.2000 

289.49 

2.000 

1 .995 

804 . 0200 

969.5000 

290.90 

3.000 

2.991 

713.9600 

871 .8000 

285.29 

3.000 

2.991 

714.0400 

073.0000 

294 . 9*3 

>♦.000 

3.900 

632.9100 

788.1000 

279.76 

4.000 

3.998 

632.9500 

787.90110 

279 . 8*3 

5.000 

4.984 

559.7000 

71 1 .8000 

273.92 

5.000 

4.904 

559.7700 

711.3000 

27i . 16 

6.000 

5.980 

493.6500 

641.7000 

268.01 

6.000 

5.980 

493.7600 

64 | .2000 

2GB. 29 

7.000 

6.976 

434. 1900 

577.5000 

261  .92 

7.000 

6.976 

434.3100 

577.9000 

2*31  .03 

8.000 

7.971 

380.7100 

519.2000 

255.46 

8.000 

7.971 

380 . 7800 

519.9000 

255.13 

9.000 

8.966 

332.6500 

466.8000 

246.26 

9.000 

8.906 

332 . C50Q 

467.3000 

247.9*3 

I0.CC0 

9.960 

289.4800 

419.1000 

240.64 

10.000 

9.9G0 

289.4300 

419.6000 

2**0 . 28 

11.000 

10.955 

250.7800 

375.4000 

232.72 

11.000 

10. 955 

250.6800 

375.9000 

232.31 

12.000 

11.949 

216.1600 

335.3000 

224. 5E 

12.000 

11.949 

216.0200 

335.7000 

224 . 1 7 

13.000 

12.942 

165. 3200 

298.3000 

216.45 

13.000 

t2.942 

185.1500 

298.5000 

216.99 

14.000 

13.936 

I57.96C0 

263.9000 

208.5! 

14.000 

13.936 

157.7700 

264 . 1 000 

200. 18 

15.000 

14.929 

133. ♦ CO 

232.0000 

201.00 

15.000 

14.929 

133.6400 

232.1000 

200.51 

16.000 

15.921 

112.7300 

202.3000 

194. 12 

16.000 

15.921 

1 12.5300 

202. 1000 

103.97 

17.000 

16.914 

94  .4710 

173.6000 

139.62 

17.000 

16.914 

94 . 3530 

172.2000 

190.00 

18.000 

17.9C6 

79.0620 

144.1000 

191 .10 

18.000 

17.906 

79.0620 

143.0000 

132. 'v> 

19.000 

18.898 

6C.4250 

1 16.9000 

198.02 

19.000 

10. 8. 10 

66.4690 

1 16.9000 

19U.il  I 

20.000 

19.899 

56. 1320 

95.6900 

204 . 35 

20 . 000 

19.809 

56. 1440 

96.2100 

203  89 

21  .000 

20.881 

47.6320 

79.7400 

208.11 

21.000 

20.881 

47.6130 

79.8900 

207.61 

22.000 

21 .eoi 

40.5340 

66.7900 

21  1 .40 

22 . 000 

21 .871 

40.5030 

65.8700 

21  1.08 

23.000 

22.852 

34.5750 

56.2200 

214.24 

23.000 

22.862 

34.5380 

56.2800 

21  7.  Oil 

24.000 

23.652 

29.5510 

47.5000 

216.74 

24.000 

23.652 

29.5120 

47.5000 

216.45 

25.000 

24.842 

25.3000 

40.2900 

218.74 

25 . 000 

24.842 

25.2630 

40.2600 

219.64 

26.000 

25.832 

21.6920 

34  2400 

220 . 70 

26.000 

25.832 

21  .6600 

3*1.1800 

280.73 

27 . 000 

26.821 

18.6230 

29.1600 

222.51 

27 . 000 

26.821 

18.5900 

29.0900 

282 .  '77 

28.000 

27.810 

16.0090 

24.8600 

224.32 

20.000 

27. BIO 

15.9930 

2*1.7500 

285.05 

29.000 

28.799 

13.7800 

21 .2300 

226.09 

29 . 000 

28.799 

13.7730 

21 .1300 

827.10 

30. COO 

29.787 

U.874g 

MARCH 

18.16)0 

227.59 

37.000 

29.787 

1  1 . 8704 

APRIL 

16.0 /00 

i 89 . 05 

TABLE 

STATION 

IV.  3 
-  912170 

HYDROSTATIC  MODEL  ATMOSPHERE, 
TAGUAC 

TABLE 

STATION 

IV.  4 
-  912170 

HYDROSTATIC  MODEL  ATMOSPHERE. 
TAGUAC 

Z 

CEO.  HT. 

P 

0 

TV 

2 

GEO.  HT. 

P 

D 

TV 

KM 

KM 

MB 

G/M3 

OEO  K 

KM 

KM 

MB 

G/M3 

DEG  K 

.000 

.000  1012.9000  1168.0000 

302.08 

.000 

.000  1012.4000  1164.0000 

302.96 

.111 

■111  1000.3000  1156.0000 

301 .38 

.111 

■  111 

999.8600  1153.0000 

302.22 

1.000 

.997 

903.5800  1069.0000 

294.48 

1.000 

.997 

903.4200  1066.0000 

295.15 

2.000 

1  .995 

004 . 0300 

968.7000 

289.14 

2.000 

1.995 

804 . 0600 

967.3000 

209.57 

3.000 

2.991 

714.0900 

872.9000 

284.99 

3.000 

2.991 

714.1300 

874.0000 

284.65 

4.000 

3.908 

632.9900 

787.9000 

279.06 

4.000 

3.988 

632.9300 

789.0000 

279.47 

5.000 

4.984 

559.6200 

71 1 .1000 

274.25 

5.  COO 

4.984 

559.6000 

71 1 .8000 

273.95 

6.000 

5.900 

493.8400 

64 I . | 000 

268.37 

6.000 

5.980 

493.6600 

641 .3000 

268. 14 

7.000 

6.976 

434.4100 

577.4000 

262.09 

7.000 

6.976 

434.1900  • 

577.7000 

261.81 

8.000 

7.971 

380.9200 

519.4000 

255.49 

8.00C 

7.971 

380.6700 

519.8000 

255.12 

9.000 

0.966 

332.8**00 

466.9000 

248.35 

9.000 

8.966 

332.5600 

467.2000 

247.99 

10.000 

9.960 

289.6700 

419.2000 

240.74 

10.000 

9.960 

289.3600 

419.3000 

240.39 

11.000 

10.955 

250.9500 

375. 5000 

232.81 

11.000 

10.955 

250.6300 

375.7000 

232.40 

12.000 

1 1 .949 

216.3300 

335.5000 

224.61 

12.000 

1 1 .949 

216.0000 

335.4000 

224.37 

17.700 

12.9*2 

185.4600 

298.5000 

216.48 

13.000 

12.942 

185. 1600 

298.1000 

216.38 

14 .000 

13.936 

158.0900 

264.2000 

208.45 

14.000 

13.936 

157.8300 

263.7000 

208.51 

15.000 

14.929 

133.94C0 

232.3000 

200 . 90 

15.000 

14.929 

133.7400 

231  .4000 

201 .31 

1C. 000 

15.921 

1 12.8000 

802.5000 

194.03 

16.000 

15.921 

1 12.7200 

201 .1000 

166.27 

17.000 

16.914 

94 . 6500 

173.3000 

190.06 

17.000 

16.914 

94.61  10 

171 . 7000 

191  .US 

IB. 000 

17.906 

79.  1910 

143.4000 

192. 35 

10.000 

17.906 

79.3300 

14  3.31)00 

198.06 

I'J.000 

18.898 

66.5000 

1 17.0000 

199.32 

19.000 

18.898 

66.7170 

1 17.1000 

198.51 

20 . 000 

!9.8ri9 

56.2530 

96.2300 

203.53 

20.000 

19.889 

56 . 3990 

95.9700 

204 . 73 

21 .000 

20.831 

47.7200 

79.97GQ 

203 . 05 

21.000 

20.891 

47.8900 

79.6000 

209 . 38 

22 . 000 

21  .871 

40,6060 

66.9900 

21  1  .29 

22.000 

■>1  .Cl 

40.7030 

66.8100 

212.72 

23 . 000 

22.862 

34.6370 

56.2700 

214.42 

23.000 

22.862 

34.0320 

56.2700 

216  66 

24 . 000 

23.652 

29.6110 

47, 67  Oil 

217.16 

24 . 000 

23.652 

29.8050 

47  .5. '00 

210.5! 

25.000 

24.842 

25.3610 

4G.6b00 

219.46 

85.000 

24.842 

25.5570 

40  2300 

261 . 31 

26.000 

25.832 

21 .7590 

34 . IbOO 

221  .93 

26 . 000 

25.832 

21 .95 CO 

34.1500 

223.95 

27.000 

26.821 

18.6990 

29.0600 

224. 15 

27.000 

26.B2I 

18.8970 

29.0500 

226 . 65 

28.000 

27.810 

16.09*0 

24.7000 

226.27 

28.000 

27.810 

16.2900 

24.8400 

2c'0.**5 

29.000 

28.799 

13.8700 

21 .2000 

2?"’ .  07 

29.000 

28.799 

14.0590 

21 .2700 

237., -*8 

30.000 

29.797 

11 .5601 

18.1400 

Pc  9 . 83 

30.000 

29.797 

12. 1502 

18.2200 

2*8. 33 

90 


MAY 


JUNE 


tadlc 

station 

/ 

KM 

.000 

.111 

1.000 

2.000 

3.000 

4.0Q0 

5 .  000 

6.  COO 
7.000 
8.000 
9.GC0 

10.000 
11.000 
12.000 
13.000 
14.000 
15.000 
16.000 
17.000 
10.000 
19.000 
20.000 
21 .COO 
22.000 
23.000 
24.000 
25.000 
25.000 
27.000 
23.000 
29.  COO 
30.003 


IV.  5 
-  912170 
OCO.  HT 
KM 


2. 991 


HY0R0STAT1C  MODEL 
TAGUAC 


ATMOSPHERE . 


7.971 


15.921 

16.91V 


19.899 
20.831 
21 .871 
22.862 
23 . 052 

24.84? 

25  832 
26.821 
27.810 
28.799 
29  7  37 


TABLE 

STATION 

Z 

KM 
.000 
.111 
1 .000 
2.000 
3.000 
4.000 
5.000 
6  000 
7.000 
8  000 
9.000 
10.000 
11.000 
u-.ooo 
13.000 
IV. 000 
15.000 
16.000 
17.000 
18.000 
19. 000 
20 . 000 
21.000 
22.000 
23.000 
2*  .  000 
25 . 000 
26  000 
27 . 000 
213.000 
.'’9.000 
30.000 


IV.  7 
*  912170 
GEO.  HT. 
KM 

.000 
.111 
.997 
1.995 
2.991 
3.988 
V  ,98V 
5.980 
6.976 
7.971 
8.966 
9.960 
10.955 
1 1 .9vg 
12.9V2 
13.935 
IV. 929 
15.921 
16.91V 
17.906 
18.898 

19.809 

20.831 
21  .671 
22.852 
23.552 
24.842 

25.832 
26.821 

27.810 
20.799 
29.707 


p 

D 

TV 

M9 

G/M3 

DEG  V. 

1  101 1 .8000 

!  161.0000 

303.52 

999.2900 

1 I49.C000 

302.85 

'  903.1400 

1063.0000 

296.0? 

>  804 . 0600 

964.9000 

290 . 30 

714  .f-)00 

373.0000 

205 . 05 

1  635.1100 

7(59.04:10 

279 , 25 

i  569.760 0 

71 3. 0000 

273.61 

1  493.6300 

642.2000 

267 . 76 

i  434.0900 

670.3000 

261 .50 

380.5300 

520.1000 

254 . 90 

.  352.4000 

4157.4000 

247.70 

1  289.1000 

419.6000 

240. 10 

i  250.4300 

375.8000 

232.  IQ 

1  216.7900 

335.6000 

224.01 

184.9300 

290.3000 

215.97 

157.5900 

263.7000 

206 .  ?,? 

133.5400 

230.9000 

201.47 

112.5900 

20”  1000 

196.02 

94 . 5590 

I  "’DO 

193.53 

79.4610 

IVl.o  00 

195. !H 

G6.SE 60 

1 16.1000 

200.93 

66.7190 

95.4600 

206. 90 

48.2300 

79.7200 

210.75 

41 . 1200 

67 . 0200 

213.75 

35.1380 

56.4S00 

215.74 

30 . 0900 

47.7500 

219.52 

26.8140 

40.5400 

221  .64 

22. 1040 

34.4600 

224.29 

19.0930 

29.4000 

226.26 

16.4570 

25.0900 

228.48 

14.2060 

21  .4600 

230 . 57 

1 2 . 27  u8 

19. 4."  90 

232. 16 

JULY 

HYDROSTATIC  MOOT'  ATMOSPHERE. 

TAGUAC 

P 

0 

TV 

MB 

G/M3 

DEG  K 

1010.6000  1160.0000 

303.63 

998.0500  II 48. 0000 

302.9! 

902.1300  1059.0000 

295.69 

803.3200 

952.8000 

290.65 

713.7000 

872.1000 

285.1! 

632.5600 

789.6000 

279 . 07 

559.2000 

7 ! 3 . 3000 

273.11 

493.0500 

642.4000 

267.39 

433.5300 

570.0000 

261 .32 

380.0100 

519.6000 

254.00 

331 .9300 

466.9000 

247.68 

288 . 7500 

419  1000 

239 . 99 

250 . 0300 

375.7000 

231 .83 

215.3600 

3.35.3900 

223.4) 

18V .4800 
157.0700 
132.9V00 
1  1 1 .9600 
94.2130 
79.51 10 
67.3230 
57 . 1 720 
40.6710 
VI .6100 
35 . V  ZOO 
30 . 3720 
26.0410 
22.3610 
19.2270 
16.5540 
14,27)0 
!.?.  3143 


230.6000 
264 .4000 

231 .eooo 

199.9000 
1C5.9000 
137.3000 
114.1000 
95.3300 
0U.OCOO 
67.6400 
57.1000 
48.3530 
41 .0700 
34.9650 
29.8100 
25.4300 
21 .7600 
16  .  •  0 


215.07 

2CG.96 

199.79 

1S5.M 

197,85 

201.69 

2C5.57 

208.93 

211.78 
214.16 
216.45 
218.83 
22C.Q8 
222 . 90 
224 . 72 

226 . 79 
220 . 24 
729 . 51 


TABLE  IV.  6  HYDROSTATIC  MODEL  ATMOSPHERE. 
STATION  -  912170  TAGUAC 

D 

G/M3 


Z 

GEO.  HT. 

i 

P 

KM 

KM 

MB 

.000 

.000 

ton 

.5000 

.111 

■  111 

998 

.9400 

1 

.000 

.997 

90? 

.9400 

2 

.000 

1 .995 

803 

.9900 

3 

.000 

2.991 

714 

.2600 

4 

.000 

3.999 

633 

.0700 

5 

.000 

4.984 

559 

.6800 

6 

■  COO 

5.980 

493 

.5000 

7 

.000 

6.976 

433 

.9400 

8 

.000 

7.971 

330 

.3900 

9 

.000 

8.966 

332 

.2400 

10 

.000 

9.960 

239 

.0200 

1! 

.000 

10.955 

250. 

.2500 

12 

.000 

1 1 .949 

215 

.5800 

13 

,000 

12.942 

164, 

.6800 

14 

.000 

13.936 

157 

.2600 

15 

.000 

14.g29 

133, 

1800 

16 

.000 

15.921 

1  12. 

.2300 

17 

000 

16.914 

94, 

3080 

'8. 

,000 

17.906 

79. 

5230 

19. 

000 

18.898 

67. 

2400 

20. 

000 

19.883 

57. 

0620 

21  . 

000 

20.881 

48. 

5640 

22. 

coo 

21  .871 

41  . 

4290 

23. 

000 

22.P-.2 

35. 

4120 

2*4  . 

000 

23.050 

30. 

3240 

25 

000 

24.842 

28. 

0130 

26. 

000 

25.832 

22. 

3490 

27. 

000 

26.B21 

19. 

2290 

26. 

000 

27.810 

16. 

6670 

29. 

000 

20.799 

!“. 

29  !0 

30. 

000 

29.787 

12. 

3-01 

1061 .OOCO 
964.1  COO 
872.9000 
789.9300 

713.6000 

642.700 0 

578.4000 
520.0000 

467.4000 

419.7000 
376.0000 
335.0000 

293.8000 
264 .2000 

231 .4000 

199.7000 

167.8000 

133.8000 

114.6 000 
95.3800 
73.9600 
67. 3  ICO 
55  8600 
40.1700 
40.0/CO 
34.7900 
29.6900 
25.3300 
2! 

I9.64U0 


AUGUST 


TV 
DEG  K 
303.86 
303. 1? 

296. 52 
290. 5C 
295.05 
279.20 
273.23 
2G7.V9 
261 .38 

254.81 
247.63 
239.91 
231 .83 

223.53 

215.32 
207 . 35 

200.53 
195.76 
196.00 
199.50 

204 . 33 
208.4? 
21  1  .53 
214.43 
216.(19 
219.32 
221 .72 
223.78 
225.70 

227.82 
229.31 
230  61 


TABLE 

STATION 

IV.  B 
-  912170 

HYDROSTATIC  MOOEL  ATMOSPHERE 
TAGUAC 

Z 

GEO.  HT. 

P 

D 

TV 

KM 

KM 

MB 

G/M3 

oea  k 

.000 

.000 

1010.3000 

1160.0000 

303.37 

.11! 

997.7500 

1 149.0000 

303.56 

1 .000 

.997 

901.8200 

1059.0000 

296.73 

2.000 

1  .995 

803.0700 

962.3000 

290.74 

3.000 

2.991 

713.4900 

871 .8000 

305. \\ 

4.000 

3.908 

632.3900 

799.1000 

279. 17 

5.000 

4.984 

559.0800 

712.7000 

273.20 

6 . 000 

5.980 

492.9900 

641 .8000 

267.57 

7.000 

6.976 

433.5200 

577.4000 

261 .57 

6 . 000 

7.971 

3U0.0600 

519.0000 

255 . 1 0 

9.000 

8.966 

332.C200 

466.4000 

247.  rJJ 

10.000 

9.960 

208.8800 

418.0000 

240 . 30 

11.000 

10.955 

25 0. 1900 

375.4000 

23?. 14 

12.000 

11.949 

215.5500 

335.7000 

2c 3. 67 

13.000 

12.942 

104.6600 

293.9000 

215.22 

14.000 

13.936 

157.2300 

264 .7000 

206.94 

15.000 

14.929 

133. 070 0 

232. 1000 

199.69 

16.000 

15.921 

1 12.0700 

2CC.0C00 

1 95 . 22 

17.000 

16.914 

94.3150 

166.0000 

197.94 

18.000 

17.905 

79.6040 

137.4000 

201.83 

19  COO 

18.898 

67.4080 

114.2000 

205 . 67 

20.000 

19.889 

57.2410 

95.6400 

200 . 72 

21.000 

20.681 

48.7230 

90. 21 70 

211 .62 

22.000 

21 .971 

4 1 . 5530 

67.650‘J 

213.97 

23.000 

22 . 862 

35.5030 

57.2200 

2 1 5 . 0Q 

24.000 

23.652 

30. 37-0 

48.60C0 

?17 . 74 

25.000 

24.8-2 

260230 

v : ,?too 

219.70 

26 . 000 

25 . 032 

22 . 3230 

*5.1X00 

r ;/ }  ,og 

27.000 

26.681 

19. I860 

c'9.6600 

?£i  80 

28.000 

27.810 

16. 5080 

25. 4 BOO 

225 . 70 

29. COO 

28 . 799 

14.2210 

21 .ecco 

227.36 

30 . 000 

29  787 

12.2634 

18.6520 

229.63 

91 


SEPTL'MBCR 


OCTOBER 


IV.  9  HYDROSTATIC  MODEL  ATMOSPHERE. 


STATION  * 

912170 

TAGUAC 

Z  GEO.  HT. 

P 

D 

TV 

KM 

KM 

MB 

C/M3 

DEO  K 

.000 

.000  1010.5000  1160.0000 

303.54 

.111 

.111 

997.9300  1148.0000 

302.04 

1 .000 

.997 

902.0400  1058.0000 

296.89 

2.000 

1 .995 

803.3300 

961 .8000 

290.96 

3.000 

2.991 

713.7800 

871 .5000 

285.33 

4. 030 

3.988 

632.6900 

789.2000 

279 . 29 

5.000 

4.984 

553.3600 

713.2000 

273.24 

6.000 

5.990 

493.2200 

642.4000 

267.45 

7.000 

6.976 

433.6900 

578.0000 

261.39 

8.000 

7.971 

380. 1600 

519.6000 

254.87 

9.000 

8.966 

332.0700 

466.9000 

247 . 76 

10.000 

9.960 

288.8900 

419.1000 

240. 14 

11.000 

10.955 

250.1800 

375.7000 

232.01 

18.000 

1 1  .949 

215.5300 

335.9000 

223.56 

13.030 

12.942 

104.6300 

298.9000 

215. 18 

14.000 

13.936 

157.2100 

264.6000 

207 . 06 

15.000 

14.929 

133.0700 

231 .8000 

200.00 

16  COO 

15.921 

1 12.0900 

199.9000 

195.36 

17.000 

16.914 

94  2810 

167.2000 

196.47 

18.000 

17.906 

79.4740 

138.2000 

200 . 30 

19.000 

18.898 

67.2210 

114.6000 

204.38 

20.000 

19.889 

57.0320 

95.6300 

207.76 

21 .000 

20.88! 

48.5140 

80.1400 

210.90 

2? . COO 

21  .871 

41 .3540 

67.6000 

213.12 

23.000 

22.962 

35.3090 

57.1700 

215. 16 

24.000 

23.652 

30.1950 

48.43G0 

217.22 

25.000 

24.842 

25 . 6630 

41  .02CC 

219.63 

26  000 

25.832 

22. 1903 

34.8600 

221 .74 

27 . 000 

26.621 

19.0670 

29.6600 

223.95 

28.000 

27.810 

15.4090 

25.2900 

226 . 07 

29.000 

28.799 

14.1423 

21 .5700 

288.2 3 

30.000 

29.707 

12.2050 

18.4900 

229.69 

NOVEMBER 

TABLE 

IV.  11 

HYDROSTATIC  mooel 

ATMOSPHERE 

STATION 

*  912170 

taguac 

Z 

GEO.  HT. 

P 

0 

TV 

KM 

KM 

MB 

GW 

DEG  K 

.000 

.000 

1010.4000 

1160.0000 

303.57 

.111 

-111 

997.8700 

1 140.0000 

302.88 

1.000 

.997 

901.9400 

1059.0000 

295.56 

2.000 

1  .995 

803. 1700 

962.0000 

290.85 

3.000 

2.991 

713.6500 

870 .9000 

295.47 

4.000 

3.988 

632.6600 

787.7000 

279.80 

5.000 

4.584 

553.4930 

71 1  .4CC0 

273.9° 

6.000 

5.980 

493.4000 

641.3000 

2C8.C6 

7.000 

6.976 

434.0300 

577.4000 

261 .89 

8.000 

7.971 

380.5600 

519.3000 

255.29 

9.000 

8.966 

332.4900 

466.7000 

248. 19 

10.000 

9.960 

289.3300 

418.9000 

240.6! 

11.000 

10.955 

250. G500 

375.3000 

232.66 

12.000 

1 1 .949 

2I6.C500 

335. 2000 

224.53 

13.000 

12.942 

185.2200 

298.1000 

216.45 

14.000 

13.936 

157.8800 

263.8000 

208.51 

15.000 

14.929 

133.7830 

231 .6000 

201 .22 

16.000 

15.921 

1 12.7200 

201 .6000 

194.79 

17.000 

16.914 

94.5630 

172.2000 

191 .36 

18.000 

17.906 

79.2830 

143.0000 

193  13 

19.000 

18.898 

66.7133 

116.5000 

199.46 

20.000 

19.889 

56.4340 

95.7200 

205.38 

21.000 

20.881 

47.9310 

79.7700 

209.31 

22.000 

21 .871 

40.8203 

66.9600 

212.31 

23.000 

22.862 

34.8430 

56.4000 

215.21 

24 . 000 

23.652 

29.8040 

47.6500 

217.87 

25 . 000 

24.842 

25.5420 

40. 3600 

220.49 

26.000 

25.832 

21.9310 

3" .2500 

223.08 

27.000 

26.821 

18.8610 

29. 1800 

225. 16 

28.000 

27.810 

16.2450 

24 ,89?0 

227.35 

29.000 

2a . 79 J 

14.0110 

PI  .2700 

229.4= 

30.000 

29.737 

12.0947 

19.061/0 

230.03 

TABLE 

STATION 

IV.  10 
-  912170 

HYDROSTATIC  MODEL 
TAGUAC 

ATMOSPHERE 

Z 

CEO.  HT. 

P 

D 

TV 

KM 

KM 

MB 

GW 

DEG  K 

.000 

.000 

1010.4000 

1159.0000 

303.61 

■  111 

■  111 

997.9100 

1 148.0000 

302 . 90 

1.000 

.997 

902.0300 

1059. 0000 

296.05 

2.000 

1  .995 

803.2900 

962.1000 

290.87 

3.000 

2.991 

713.7400 

871 .4000 

285.34 

4.000 

3.986 

632.6900 

738. 6000 

279.49 

5.000 

4.984 

559.4200 

712.4000 

273.57 

6.000 

5.980 

493.3500 

641.6000 

267.86 

7.000 

6.976 

433.8800 

577.6000 

261 .67 

8.000 

7.971 

380.3800 

519.5000 

255.08 

9.000 

8.966 

332.2900 

466.9000 

247.93 

to. 000 

9.960 

289.1100 

419.3000 

£40.21 

11.000 

10.955 

250.3800 

375.7030 

232.19 

12.000 

11 .949 

215.7400 

335.7000 

223 . 87 

13.000 

12.942 

104.8600 

298.6000 

215.63 

14.000 

13.936 

157.4700 

264.2000 

£07.62 

15.000 

14.929 

133.3600 

231  .4000 

200 . 76 

16.000 

15.921 

112.3700 

200.5000 

195.22 

17.000 

16.914 

94.3810 

169.9000 

193.5! 

18.000 

17.906 

79.3340 

140.4000 

196.81 

19.030 

18.699 

66.9490 

115.2000 

202.4? 

20.000 

I9.BB9 

56.7280 

95.6900 

206.53 

21 .000 

20.881 

48.2190 

79.9200 

210.19 

22.000 

21 .671 

4 1 . 0930 

67. 1603 

£13. !5 

23.000 

22.862 

35.0970 

56.S0G0 

£15. 16 

24.000 

23.652 

30.0060 

48.0GU0 

£17.53 

25.000 

24.842 

25.7100 

40.7000 

££0 . 07 

26.000 

25.832 

22.0700 

34.5100 

£££.00 

27 . 000 

26.821 

18.9780 

29.3700 

££5. I£ 

20.000 

27.810 

16.3440 

25.07C0 

££7. 13 

29.000 

28.799 

14.0900 

21  .4800 

£29 . 21 

30.000 

29.787 

!2. 1731 

DECEMBER 

18.3300 

231 .29 

TABLE 

STATION 

IV.  12 
•  912170 

HYDROSTATIC  MODEL 
TAGUAC 

ATMOSPHERE 

Z 

GEO.  HT. 

P 

D 

TV 

KM 

KM 

MB 

GW 

DEG  K 

.000 

.000 

1012.0000 

1164.0000 

302.85 

.111 

.111 

999.4300 

1152.0000 

302. 18 

1.000 

.997 

903. 1000 

1064.0000 

295.65 

2.000 

1  .995 

803.9600 

965.0000 

290.22 

3.000 

2.991 

714 .24*00 

871 .9000 

285.38 

4.000 

3.988 

633.2000 

788.0000 

279.93 

5.000 

4.984 

559.9800 

71 1 .9000 

274 . 02 

5.000 

5.980 

493.9200 

641 .9000 

268.07 

7.000 

6.976 

434.4100 

578.2000 

261.75 

8.000 

7.971 

380.6600 

519.9000 

255. 19 

9.000 

8.966 

332.7300 

467.4000 

248.00 

10.000 

9.960 

289.5100 

419.7000 

240.33 

11.000 

10.955 

250  7500 

375.9000 

232 . 40 

12.000 

1 1  .949 

216. 1000 

335.7000 

224.2? 

13.000 

-.2.942 

185.2200 

298 .6000 

216. 12 

14.000 

13.93b 

157.8400 

264.1000 

208.20 

15.000 

14 .929 

133.7100 

£32.0000 

200.81 

16.000 

15.921 

1 12.6100 

802. 1000 

194 . 14 

17.000 

16.914 

94 . 3770 

I73.3JC0 

109.72 

16.000 

17.906 

78.9980 

1W8000 

191  37 

19.000 

18.898 

66.3910 

1 16.6000 

198.41 

20.000 

19.889 

56. 1220 

95.4900 

204 . 75 

21.000 

20.88! 

47.8480 

79.4800 

209 . 00 

22.000 

2! .871 

40.5760 

66.5500 

212. *9 

23.000 

22.862 

34.6410 

55.9700 

215.59 

24.000 

23.652 

29.6370 

47.3500 

?  1 8 .  OG 

25.000 

24.842 

25 . 3990 

40.1700 

820 . 26 

26 . 000 

25.832 

21 .8000 

34.1600 

222.30 

27 . 000 

26.821 

18.7390 

29.0810 

224 . 47 

28.000 

27.810 

16.1310 

24.8?u0 

226 . 39 

89 . 000 

28.799 

13.9040 

21  2100 

828  Vi 

30 . 000 

29  787 

1  1  .64991 

'0.  >'•  M 

2=9.  .'5 

92 


AM*UAL 


TA9LE 

STATION 

Z 

KM 

.000 

•  III 
1.000 
2.000 
3.000 
* .  000 
5.000 
6.000 
7.000 
8.000 
9.000 
10.000 
11.000 
12.000 
1 3 .  COO 
l*.0C0 
15.000 
16.000 
17.000 
18.000 
19.000 
20.000 
21  .000 
22.000 
23.000 

2*i .  000 

25 . 000 
26.000 
27.000 
29.000 
29-000 
30.000 


tV.  13  HYDROSTATIC  MODEL  ATMOSPHERE. 
«  912170  TAGUAC 


CEO.  HT.  P 


KM 

MB 

.000 

1011.5000 

-III 

998.9600 

.997 

902.7*00 

1.995 

803.6700 

2.991 

713.9*00 

3.988 

632.8*00 

4 . 99* 

559.5800 

5.980 

*93.5000 

6.976 

*3*. 0100 

7.971 

380.5000 

8 .956 

332. *000 

9. 990 

289.2100 

10.955 

250. *300 

11. 9*9 

215.0500 

12.9*2 

16*. 9900 

13.935 

157.6000 

1* .929 

133. *600 

15.921 

1 12. *500 

16.01* 

9*. *2*0 

17.906 

79.3210 

18.898 

66.B700 

19.889 

56.6230 

20.881 

*8.1150 

21 .871 

*0.9920 

22.PG2 

3*. 9990 

23.o52 

29.9390 

2*. 8*2 

25  6560 

25.832 

22.0230 

26.821 

18.9330 

27.810 

16.3000 

23.799 

l* .0510 

25.787 

12.1277 

D 

TV 

G/M3 

DEO  K 

1163.0000 

303.05 

1 151 .0000 

302.3* 

1063.0000 

295.68 

96*. 8000 

290.19 

872.2000 

285. 1* 

788.7000 

279.51 

712.3000 

273. G9 

6*1 .8000 

267.87 

577.8000 

201 .67 

519.6000 

255.10 

*67.0000 

2*7.90 

*19.3000 

2*0 . 20 

375.7000 

232.2 9 

325.6000 

22* . u* 

298.5000 

215.85 

cl*. 1COO 

20? . 66 

231  .7000 

200.6? 

201 .0020 

IS-*.  93 

!70  2000 

193.23 

1*1  .2000 

195.70 

115.8000 

201 . 17 

95.7000 

206.13 

79.9000 

209.30 

67.1200 

212.77 

56.6000 

215. *3 

*7.8700 

217. P8 

*0.5000 

220.25 

3*. *800 

222.52 

29.3700 

22*. 60 

25.0500 

226.66 

21  .*200 

228.50 

18.3600 

230.08 
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TABLE 

IV.  13 

HYDROSTATIC  MOOCL  ATMOSPHERE 

STATION 

-  912170 

TAGUAC 

Z 

GEO.  HT. 

P 

0 

TV 

KM 

KM 

MB 

G/M3 

DEO  K 

.000 

.000 

1011.5000 

t 163.0000 

303.05 

•  111 

•  111 

998.9600 

1 151 .0000 

302.34 

1.000 

.997 

902.7400 

1063.0000 

295.88 

2.000 

1.995 

803.6700 

964.8000 

290. 19 

3.000 

2.991 

713.9*00 

872.2000 

285.14 

<*.000 

3.988 

632.8400 

788.7000 

279.51 

5.000 

>♦.994 

559.5800 

712.3000 

273. G9 

6.000 

5.980 

493.5000 

641.8000 

267.87 

7.000 

6.976 

434.0100 

577.8000 

261 .67 

B.000 

7.971 

380.5000 

519.6000 

255. 10 

9.000 

8.906 

332.4000 

467.0000 

247.99 

10.000 

9.950 

289.2100 

419.3000 

240. 30 

1 1 .000 

10.955 

250 . 4900 

375.7000 

232.29 

l?.000 

11.949 

215.8500 

335.6000 

224 . 04 

it.  coo 

12.942 

164.9800 

298.5000 

215.85 

I4.0C0 

13.935 

157.6000 

204 . 1000 

2C7.  &5 

15.000 

14.929 

133.4800 

231 .7000 

eca .67 

16.000 

15.921 

112.4500 

201 .0020 

ISh.93 

17.000 

16.014 

94.4240 

170 .2000 

193.23 

18.000 

17.906 

79.3210 

141.2000 

195.70 

19.000 

18.896 

66.6700 

115.8000 

201 . 17 

c 0.000 

19.889 

56.6230 

95.7000 

206. 13 

21 .000 

20.881 

48.1150 

79.9000 

209.80 

22.000 

21.871 

40.9920 

67.1200 

212.77 

23.000 

22.862 

34.9990 

56.6000 

215.43 

24.000 

23 .b52 

29.9390 

47.8700 

2i7. en 

25 . 000 

24.842 

25  6560 

40.5800 

2c?0 .25 

26.000 

25.832 

22.0230 

34.4800 

222.52 

27.000 

26.821 

18.9330 

29.3700 

224.60 

28.000 

27.810 

16.3000 

25.0500 

226.66 

23.000 

28.799 

14.0510 

21.4200 

228.50 

30.000 

£8.787 

12.1277 

18.3600 

230.08 

APPENDIX  A 


EXAMPLES  OF  WIND  STATISTICS  FOR  EGLIN  AFB,  FLORIDA 


Appendix  A  gives  some  examples  of  further  computations  and  graphical  dis¬ 
plays  of  wind  statistics  that  can  be  derived  from  the  statistical  parameters 
presented  in  table  I.  These  illustrations  should  aid  the  user  of  the  RRA  to 
understand  the  functional  relationships  of  the  probability  wind  models  and, 
thus,  to  develop  an  appreciation  of  the  powerful  properties  of  the  bivariate 
normal  probability  distribution  function  (PDF). 

All  illustrations  for  this  appendix  are  derived  from  the  five  wind  compo¬ 
nent  statistical  parameters  from  table  1.1  for  January  and  table  1.7  for  July 
for  nine  selected  altitudes.  These  selected  altitudes  are  2,  4,  8,  12,  16,  20, 
24,  28,  and  30  km. 

1 .  Windspeed  (Tables  A-1  and  A-2) 

The  five  wind  components  from  table  I  are  used  as  inputs  to  the  generalized 
Rayleigh  PDF,  equation  (29),  and  numerically  integrated  as  indicated  by  equation 
(30)  to  obtain  the  PDF  for  windspeed.  The  PDF  is  then  interpolated  to  obtain  the 
percentile  values  for  windspeed,  as  shown  in  tables  A-1  and  A-2. 

2.  Frequency  of  Wind  Direction  (Figures  A-1  through  A-18) 

The  derived  frequencies  for  wind  direction  shown  in  figures  A-1  through 
A-18  were  obtained  using  the  five  wind  component  parameters  from  tables  1.1 
and  1.7  as  input  values  in  equation  (35).  The  limits  of  integration  (performed 
numerically)  are  over  the  22.5-degree  interval  for  each  of  the  16  compass  points. 
These  graphs  give  the  percentage  frequency  that  the  wind  will  blow  from  the 
direction  intervals. 

3.  Mean  Wind  Components  and  80th  Interpercentile  Range  of  Wind  Components 

(Figures  A-1 9  through  A-36)  3  C 

The  wind  component  means  with  respect  to  any  orthogonal  axes  are  obtained 
by  using  the  zonal  and  meridional  mean  wind  components  in  equations  (44)  and 
(45).  These  component  means  form  the  circle  shown  in  figures  A-1 9  through 
A-36.  Further,  the  zonal  and  meridional  wind  component  variances  and  corre¬ 
lation  coefficients  are  used  in  equations  (46)  and  (47)  to  obtain  the  variances 
with  respect  to  any  orthogonal  axes.  These  rotated  component  variances  and  the 
rotated  component  means  are  used  in  equation  (8)  to  obtain  the  80th  interper¬ 
centile  range  of  wind  components  and  are  then  illustrated  in  figures  A-1 9 
through  A-36. 

4 .  Probability  Ellipses  (Figures  A- 37  through  A-54) 

Using  the  five  wind  component  parameters  from  tables  1.1  and  1.7  and 
p  =  0.50,  p  =  0.95,  and  p  =  0.99  as  input  values  to  equation  (13),  the  wind 


probability  ellipses  shown  in  figures  A- 37  through  A-54  were  obtained  by  com¬ 
puter  graphics.  The  statistical  inferences  are,  for  example,  that  50  percent 
of  the  wind  vectors  lie  within  the  smaller  ellipse  and  99  percent  of  the  wind 
vectors  lie  within  the  outer  ellipse.  These  probability  ellipses  are  illus¬ 
trated  using  the  standard  meteorological  coordinate  system  explained  in 
section  I . B . 1 . 

5.  Conditional  Windspeed  Given  the  Wind  Direction  (Figures  A-55  through  A-72) 

The  five  wind  component  parameters  from  table  1.1  and  table  1.7  are  used 
to  evaluate  the  conditional  PDF,  equation  (41).  Interpolations  of  the  condi¬ 
tional  function  are  made  to  obtain  the  5th,  15th,  50th  (median),  85th,  95th, 
and  99th  conditional  percentile  values  of  windspeed  given  the  wind  directions 
are  as  shown  in  figures  A-55  through  A-72.  The  conditional  mean  windspeed, 
given  the  wind  direction,  is  obtained  from  equation  (40).  The  conditional 
mode  (most  probable)  windspeed,  given  the  wind  direction,  is  obtained  from 
equation  (38).  The  conditional  mean  windspeed  and  the  conditional  windspeed 
modal  value,  given  the  wind  direction,  are  also  shown  in  these  figures.  For 
some  figures,  the  conditional  windspeed  values  are  invalid  for  the  given  wind 
direction  near  270°  (from  the  west).  This  is  caused  by  the  lack  of  computa¬ 
tional  precision  in  evaluating  equations  (40)  and  (41)  when  the  arguments 
for  the  Gaussian  probability  distribution  have  large  negative  values,  i.e., 
when  the  coefficients  (b/a)  become  less  than  -4  in  these  equations. 

This  appendix  contains  only  a  few  of  the  many  options  in  presenting  wind 
statistics  illustrations. 


TABLE  A-l.  DERIVED  (RAYLEIGH)  PERCENTILES  FOR  WINDSPEED, 
JANUARY,  TAQUAC,  GUAM  ISLAND 


Altitude  (km 
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JULY,  TAQUAC,  GUAM  ISLAND 


Altitude  (km 
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APPENDIX  B 


RANGE  SPECIFIC  INFORMATION  AND  THERMODYNAMIC  QUANTITIES  FOR 

TAQUAC,  GUAM  ISLAND 

1.  Range  Specific  Information 

To  prevent  further  character  size  reduction  for  tables  I  through  IV,  cer¬ 
tain  range-specific  information  has  been  omitted.  This  important  information 
is  given  in  tabl e  B-l . 

TABLE  B-l 

Header  Record  0-30  Km 


Table  Number - 0 

Data  Source  (1  =  DATSAV,  2  =  WDC-A) - 1 

Call  Letters - PGAC 

WMO  Number . - . 912170 

Latitude - 13.33 

Direction  (N  or  S) - N 

Longitude - - - 144.51 

Direction  (E  or  W) - E 

Elevation  in  Meters - 111 

Start  Period  of  Record  (Mo-Yr ) - 173 

End  Period  of  Record  (Mo-Yr) - 1279 

No.  of  Time  Windows  (0,  1  or  2) - 2 

Start  Time  Window  #1  (Hr-MNZ) . 2200 

End  Time  Window  #1 - 200 

Start  Time  Window  #2 - 1000 

End  Time  Window  #2 - 1400 

Date  of  RRA - 1080 

Altitude  Range  of  RRA  Low  Level  (Km) - 0 

Altitude  Range  of  RRA  High  Level  (Km) - 30 

Standard  Deviation  of  Thermodynamic  Limits - - --6.0 

Wind  Limits - 6.0 


2.  Thermodynamic  Quantities 

This  section  presents  examples  of  further  computations  and  graphical  dis¬ 
plays  of  pressure,  density,  and  virtual  temperature  statistics  that  can  be 
derived  from  data  given  in  tables  II,  III,  and  IV.  No  attempt  is  made  to 
present  complete  nor  exhaustive  illustrations  that  can  be  made  to  aid  in  vis¬ 
ualizing  the  relationships  that  can  be  made  from  the  data  in  tables  II  and  IV. 
The  choices  are  those  that  aided  the  committee  to  verify  the  reasonableness 
of  the  tabulations. 

2.1  Monthly  Means  from  the  Annual  Mean 

The  hydrostatic  model  values  in  table  IV  are  used  to  compute  (1)  the 
monthly  mean  differences  relative  to  the  annual  mean  values  of  pressure. 
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2 .  i  Monthly  Means  from  the  Annual  Mean 

The  hydrostatic  model  values  in  table  IV  are  used  to  compute  (1)  the 
monthly  mean  differences  relative  to  the  annual  mean  values  of  pressure,  den¬ 
sity,  and  virtual  temperature  expressed  in  percent  and  (?)  the  monthly  mean 
difierence  in  virtual  temperature  for  the  annual  mean  virtual  temperature 
expressed  in  decrees  Kelvin.  Examples  of  these  four  statistics  are  given  in 
table  B-2  for  January  and  table  B-3  for  July.  Graphical  displays  of  the  four 
statistics  contained  in  tables  B-2  and  B-3  are  shown  in  figures  R-i  through 
B-8.  Also,  the  relative  differences  between  the  monthly  mean  values  from 
table  IV-1  through  IV-12  for  all  months  from  the  annual  mean  values  (table 
IV- 1 3 )  are  illustrated  in  figure  B-9  for  pressure,  in  figure  B-10  for  density, 
and  in  figure  B-ll  for  virtual  temperature.  The  monthly  mean  virtual  tempera¬ 
ture  differences  from  the  annual  mean  virtual  temperature  for  all  months  are 
given  in  figure  B-12.  The  simple  sum  of  the  monthly  mean  differences  from 
the  annual  mean  values  of  these  quantities  is  not  zero.  This  is  because  the 
annual  mean  statistical  parameters  are  computed  (see  section  III.  C.3  of  text) 
by  weighting  the  monthly  means  by  the  number  of  observations  of  each  month. 

2 . 2  Coefficients  of  Variation  and  Derived  Correlation  Coefficients 

The  coefficient  of  variation,  Cy,  is  defined  by  the  standard  deviation 

with  respect  to  the  mean  divided  by  the  mean.  The  coefficients  of  variation 
for  pressure,  CyP,  and  density,  CyD,  were  computed  using  the  standard  devia¬ 
tions  from  table  II  and  the  hydrostatic  mean  values  from  table  IV.  The  coeffi¬ 
cient  of  variation  for  temperature  uses  the  standard  deviations  of  virtual 
temperature  from  table  III  to  the  altitude  where  virtual  temperature  exists. 
Above  this  altitude,  the  standard  deviations  of  temperature  are  from  table  II. 
The  mean  values  for  temperature  (virtual  temperature  to  the  altitude  where  it 
exists)  are  taken  from  table  IV.  No  distinction  is  made  in  the  table  headings 
in  table  B-4  (January)  and  table  B-5  (July)  and  all  related  figures  between 
virtual  temperature  and  temperature. 

from  the  coefficients  of  variation  for  pressure,  density,  and  temperature 
(virtual  temperature  to  the  altitude  where  it  exists),  the  correlation  coeffi¬ 
cients  between  these  quantities  are  derived  using  Buell's  method  (see  reference 
in  text).  The  equations  for  these  derived  correlation  coefficients  are 


KIM) 


(CVT)2  +  (C  VP):  -  (CVD)2 

2icyr  typ| 


( B-l 


r(  I’.f )) 


(CVD)2  -  (CVT)2  +  (CVP)2 
2[CVD  •  C  vPi 


(B-2 


-  <('VD)-  -  (CVT) 
ICVT  (' vuJ 


The  correlation  coefficients  in  tables  13-4  and  B-5  are  derived  from  the 
above  equations. 

A  test  for  the  validity  of  the  derived  correlation  coefficients  is  that 
all  three  of  the  following  inequalities  be  satisfied. 

cvp-Icvd  +  cyn  <0 

CVD  -  [CVT  +  CyP]  <  0 

cyr  -  [cvp  +  cvDi  <n 


In  these  examples  (tables  B-4  and  B-5)  the  numerical  values  from  equation  (B-4) 
are  all  negative;  hence,  the  derived  correlation  test  is  considered  valid.  The 
rare  exceptions  to  this  test  for  several  RRAs  occur  at  the  extreme  highest  alti 
tudes,  where  samples  sizes  for  the  statistical  sample  are  small. 

The  statistical  parameters  from  table  B-4  (January)  and  table  B-5  (July) 
are  illustrated  in  figures  B-13  through  B-16. 

For  all  months  the  CyP  values  are  shown  in  figure  B-17,  the  CyD  values 

are  shown  in  figure  B-18,  and  CyT  values  are  shown  in  figure  B-19.  If  the 

abscissa  on  the  figures  for  the  coefficient  of  variation  were  multiplied  by 
100,  these  figures  would  show  the  percentage  of  the  random  dispersion  of 
these  quantities  over  the  month  with  respect  to  the  monthly  mean  for  these 
thermodynamic  quantities. 

The  derived  correlation  coefficients  for  all  months  are  illustrated  in 
the  following  figures: 

a)  Figure  B-?0  gives  r(P,D). 

b)  Figure  B-21  gives  r(P,T). 

c)  Figure  B -22  gives  r(T,D). 
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.0152 

-.1039 

.5913 

-.0635 

-  .00-4 

-.rcec 

-.0J.r 

19.000 

.0091 

.0173 

.0149 

.0155 

.5105 

-.8519 

-.003! 

.0  360 

-.0115 

20.000 

.  0098 

.0166 

.0133 

-.0594 

.6018 

-.0330 

- .  C?C7 

- .  (106.9 

-.0105 

21.000 

.0090 

.0136 

.0115 

.1452 

.5442 

-.7530 

-.016! 

.0.170 

-.0111 

22.000 

.0095 

.0189 

.01  17 

.2731 

.  4858 

- . 7074 

-.016! 

-  .0083 

- .Oi  ;g 

23.000 

.0101 

.0124 

.0119 

.3600 

.4323 

-.6534 

-.O>0 

- . 0027 

-.0106 

24 . 000 

.0104 

.0178 

.0108 

.4861 

.4761 

-.537! 

-.01  18 

-.0105 

-.0103 

25.000 

.01  14 

.0109 

.0103 

.5049 

.5721 

-.4191 

-.GC?3 

-.0109 

-.01 00 

26.000 

.0123 

.0125 

.0097 

.3776 

.6919 

-.4073 

- .  00'i3 

-.0096 

-.0151 

27.000 

.0134 

.0133 

.0094 

.  3570 

.  7496 

-.3405 

- .  C(/J3 

-.0056 

-.0170 

29.000 

.0136 

.0135 

.  0091 

.3520 

.7745 

-.3195 

-.0099 

-.0093 

-  .  M !  B9 

29.000 

.0144 

.0148 

.01  10 

.  34  1  i 

.7162 

-.41 15 

-.0115 

-.0106 

-.oi*v 

30.000 

.0146 

.0153 

.0107 

.2903 

.7442 

-.4153 

-.01  m 

-.0100 

- .  0  ! 

TABLE  B-5 
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LEVEL 

CVP 

CVO 

CVT 

R1P.TI 

R(P,D> 

RIT.D1 

DC  VP 

DCVO 

DCVT 

.000 

.0018 

.0066 

.0064 

.0057 

.2680 

-.9619 

-.01  12 

-.0016 

-.0021 

.111 

.0018 

.0062 

.0060 

.0505 

.2432 

-.9565 

-.0104 

-.0016 

-.0020 

1.000 

.0019 

.0031 

.0025 

.0384 

.581 1 

-.7909 

- . 0C37 

-.0013 

.0025 

2.000 

.0019 

.0034 

.0030 

.1029 

.  4554 

-.8325 

-.0045 

-.0015 

—  .  0 Cl' j 

3.000 

.0019 

.0033 

.0031 

.2225 

.3783 

-.8183 

-.0044 

-.0018 

“.000! 

4.000 

.0022 

.0032 

.0030 

.2498 

.4498 

-.7525 

-.0040 

- . CO! 9 

-.000- 

5. COO 

.0023 

.0035 

.0031 

.1827 

.4830 

- . 7726 

-.0044 

-.0019 

‘.0007 

6.000 

.0022 

.0038 

.0037 

.2190 

.3683 

-.8265 

-.0053 

-.0020 

- . 

7.000 

.0025 

.0041 

.  0J43 

.3578 

.2247 

-.8295 

- . 0C59 

- . 0026 

- .  00?  5 

8.000 

.0027 

.0042 

.  0047 

.4594 

.  1229 

-.8251 

- . CC62 

- .003! 

- . 000? 

9.000 

.0030 

.0039 

.0049 

.5983 

.0175 

-.7907 

- . 0059 

-.0040 

-.0000 

1 7 . 330 

.  OC  -1 1 

«r*TO 

nnav 

cn 

- r  o4un 

_  764Q 

-  nnso 

-  rv.49 

-  on?n 

11.000 

.0039 

.0037 

.0056 

.7556 

-. 101 ! 

- . 7280 

- . 0054 

-  .  CO  GO 

-  .001^ 

12.000 

.0046 

.0034 

.0050 

.8080 

-.0339 

-.6162 

- .0047 

-  .  0070 

.0000 

1  5.000 

.0055 

.0042 

.0064 

.7572 

•  Min 

-.5399 

- . 0062 

- .  onvr, 

* .00 

14 .030 

.  0060 

.  0058 

.0078 

■  6’*90 

.  1  301 

-.6396 

-.0075 

- .  ugg; 

‘  .  OG^O 

15.030 

.0067 

.0085 

.009! 

.44  5? 

.  3132 

-.7110 

-.0110 

-  .CjVJ 

OOG0 

16.000 

.0073 

.0143 

.0131 

.  1  0P4 

.  4 1 30 

-.8606 

- . 020 1 

-.0061 

-  .1  JvJb*. 

17.000 

.0074 

.0159 

.0131 

-.1386 

.  5772 

- .  B097 

-.0217 

-A)',  tb 

- .  0 !  (,t’ 

18.003 

.0076 

.0121 

.0092 

-.0277 

.6502 

-.7775 

-.01 37 

r-  v 7 

- .  0 1  CG 

19.030 

.0071 

.0104 

.0008 

.1512 

.6620 

-  .74  08 

-.0121 

0  .k  G 

-.ors; 

2O.000 

.0078 

.0101 

.  3090 

.2862 

-  6720 

-.01 13 

- .  r ,  ■  ■  *  <  * 

21.000 

.0077 

.  0097 

.  11078 

.2218 

.  h :  it 

-  .6338 

- . 0096 

0 

*.0M  !G 

22.000 

.0080 

.01  1  1 

.  0090 

.1493 

A  /V* 

-  .  7C20 

-.0121 

- .  vyj 

-  .1101 

23.000 

.  0083 

•  Cl  15 

.0099 

.  1269 

.c -r* 

-  6904 

-.Cl  19 

-  v.* 

- .  J !  1  ! 

24 . 000 

.0090 

•  Oil  1 

.0096 

.2042 

.  6G  1  b 

-.6095 

- .  c :  C3 

•  Wf  ‘G 

-  .  0  !  !  G 

25.000 

.0090 

.01  I  1 

.0052 

.  1735 

.fO!Q 

- .  r  ^  3! 

-.O'  .0! 

.  o )  1  - 

26.000 

.0095 

.0!  14 

.  0099 

.2278 

.evo 

-.6j:9 

- . p: -y 

.:  *r 

-  -  01*  -T 

27.000 

.0100 

.0105 

.0098 

3762 

.  *  '  ?0 

-.4-66 

-.0-73 

-.on; 

28.000 

.0106 

.0111 

.0796 

.  3990 

6'  U* 

-.431! 

-  .01-21 

-  .0.'  0. 

. -'ll.  ) 

29.000 

.0117 

.01 1 1 

.6108 

.5174 

.  yoo 

-,4  3fi 

-.nr, ? 

-  Mis 

.V  If 

30.000 

.0127 

.0111 

.0117 

.5901 

.  G,  Vi 

- .  3797 

-.0:  * i 

.0!  .<i 

.0:01 
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